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EXISTENCE OF GLOBAL STRONG SOLUTIONS TO A BEAM-FLUID 

INTERACTION SYSTEM 

C. GRANDMONT* * AND M. HILLAIRETt 


Abstract. We study an unsteady non linear fluid-structure interaction problem which is a 
simplified model to describe blood flow through viscoleastic arteries. We consider a Newtonian 
incompressible two-dimensional flow described by the Navier-Stokes equations set in an unknown 
domain depending on the displacement of a structure, which itself satisfies a linear viscoelastic 
beam equation. The fluid and the structure are fully coupled via interface conditions prescribing 
the continuity of the velocities at the fluid-structure interface and the action-reaction principle. 
We prove that strong solutions to this problem are global-in-time. We obtain in particular that 
contact between the viscoleastic wall and the bottom of the fluid cavity does not occur in finite 
time. To our knowledge, this is the first occurrence of a no-contact result, but also of existence 
of strong solutions globally in time, in the frame of interactions between a viscous fluid and a 
deformable structure. 


1. Introduction 

In this paper, we focus on the interactions between a viscous incompressible Newtonian fluid 
and a moving viscoelastic structure located on one part of the fluid domain boundary. This work 
is motivated by the study of blood flow in arteries and the fluid-structure interaction (FSI) model 
we consider here can be viewed as a simplified version of a standard model/benchmark for FSI 
problems/solvers in hemodynamics [421 . [131 . Here we are interested in global existence result 
and the possibility of collapse of the arterial wall. Consequenlty, we investigate whether or not 
collision occurs between the moving boundary and the bottom of the fluid cavity together with 
the existence of global-in-time strong solutions. 

A vast majority of studies on the existence of solution for fluid-structure interaction problems 
concerns the motion of a rigid solid in a viscous incompressible Newtonian fluid, whose behavior is 
then described by the Navier-Stokes equations (see [J5J [351 13H [HI HI [13 HH11321 EH HU EE1 HZ] ) ■ A 
challenging point is the possibility of body-body or body-boundary collisions. In particular, these 
existence results are valid up to contact, except those of [44] and jT4j where special weak solutions 
after possible collisions are built in the 2D case and in the 3D case respectively. The contact 
issue is further investigated in [28] and [29] where a no-collision result is proven in a bounded 
two-dimensional cavity. The three-dimensional situation is also explored in [ 3T| and in [301 . We 
mention that, in [46l 144] . collisions, if any, are proved to occur with zero relative translational 
velocity as soon as the boundaries of the rigid bodies are smooth enough and the gradient of the 
surrounding fluid velocity field is square integrable. A complementary study of the influence of the 
smoothness of boundaries on the existence of collisions has been recently tackled in [17] and [18]. 
In [12] it is proved that below a critical regularity of a 2D rigid body boundary, namely C 1,1 / 2 , 
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collisions may occur ; in m and [20] slip boundary conditions at the fluid-structure interface are 
introduced and existence of weak solutions up to collision/existence of finite time collisions are 
proven respectively. 

Fewer studies consider the case of an clastic structure evolving in a viscous incompressible 
Newtonian flow. We refer the reader to [12] and [4] where the structure is described by a finite 
number of eigenmodes or to [3] for an artificially damped elastic structure while, for the case of a 
three-dimensional elastic structure interacting with a three-dimensional fluid, we mention mm] 
in the steady state case and El 0 m E3! for the full unsteady case. In the latter, the authors 
consider the existence of strong solutions for small enough data locally in time. In this field, the 
question of selfcontact and/or body-boundary collisions remains entirely open to our knowledge. 

Concerning the fluid-beam, or more generally fluid-shell coupled systems, that we consider 
herein, the 2D/1D steady state case is considered in [23 for homogeneous Dirichlet boundary 
conditions on the fluid boundaries (that are not the fluid-structure interface). Existence of a unique 
strong enough solution is obtained for small enough applied forces. In the unsteady framework 
we refer to [3 where a 3D/2D fluid-plate coupled system is studied and where the structure is a 
damped plate in flexion. The case of an undamped plate is studied in [ 23] - The previous results 
deal with the existence of weak solutions, i.e. in the energy spaces, and rely on the only tranversal 
motion of the elastic beam that enables to circumvent the lack of regularity of the fluid domain 
boundary (that is not even Lipschitz). These results also apply to a 2D/1D fluid-shell coupled 
problem which is recently considered in j4T]- I 11 this reference, the authors give an alternative 
proof of existence of weak solutions based on ideas coming from numerical schemes [21] . The 
existence of strong solutions for 3D/2D, or 2D/1D coupled problem involving a damped elastic 
structure is studied in m [39] [38]. The proofs of [39] SB] are based on a splitting strategy for the 
Stokes system and on an implicit treatment of the so called fluid added mass effect. Finally, the 
coupling of a 3D Newtonian fluid and a linearly elastic Koiter shell is recently studied in [37]. In 
this study, the mid-surface of the structure is not flat anymore and existence of weak solutions is 
obtained. 

In all the results we mentioned up to now, the existence of strong solutions is obtained locally 
in time while existence of weak solutions is obtained up to contact between the elastic structure 
and the fluid cavity. Consequently, if one wants to prove existence of solutions globally in time, 
either one should give a sense to solutions in case of collision or one should prove that no collision 
occurs. In this paper we investigate these issues on a 2D/1D fluid-beam coupled problem in 
which the beam is viscoelastic and moves only in the tranversal direction relatively to its flat mid¬ 
surface. One already knows that a unique strong solution exists locally in time (39) . whereas weak 
solutions exist as long as the beam does not touch the bottom of the fluid cavity (see [3 or [23] 
and [41j in the undamped beam case). Note that, in this model, the fluid domain is a subgraph 
whose regularity depends on the structure displacement. In connexion with the rigid-body case, 
we mention that energy estimates ensure that the beam displacement belongs to L/°(.ff 2 ) that 
embeds in L/°(C^’ Q ), for a = 1/2. This corresponds precisely to the threshold exhibited in [17]. In 
this paper, the strategy we develop is first to prove that no contact occurs and, next to propagate 
the solution regularity. In the second step, the cornerstone is an elliptic regularity result for the 
inhomogeneous Stokes system valid for nonstandard regularity of the domain boundary. 

2. General setting, main result and formal argument 

We consider a 2D container whose boundary is made of a ID viscoelastic beam, which is a 
simplified linear viscolclastic Koiter shell model m- Due to the complexity of the underlying 
fluid-structure interaction problem we assume that only the upper part of the fluid cavity is 
moving. The fluid domain denoted F(t) C R 2 depends then on time since it depends on the 
structure displacement i). It reads 

T(t) := {(x, y) e R 2 , it (0, L) , y € (0, h{x, t))} . 

where ( x , t) H>■ h(x , t) = 1 +rj(x, t ) stands for the “deformation” of the beam. We assume that the 
fluid is two dimensional, homogeneous, viscous, incompressible and Newtonian. Its velocity-field 
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u and internal pressure p satisfy the incompressible Navier-Stokes equations in F(t)\ 

(1) pf(d t u + u • Vzi) — div a(u,p) = 0, 

( 2 ) div u = 0 . 

The fluid stress tensor cr(u,p) is given by the Newton law: 

a(u,p) = p(Vu + Vu T ) — pi 2 • 

Here p denotes the viscosity of the fluid and pf its density. The structure motion is given by a 
linear damped beam equation in flexion: 

(3) p s d tt h-(3d xx h + a d xxxx h - 7 d xxt h = <j>(u,p,h), on (0 ,L), 

where a, fi, 7 are non-negative given constants and p s > 0 denotes the structure density. 

We emphasize that the beam equation is set in a reference configuration whereas the fluid 
equations are written in Eulerian coordinates and consequently in an unknown domain. Note 
moreover that we choose to write the beam equation on h and not on the beam displacement 
77 = h — 1 as is standard, since this is equivalent in the case one considers here and since it 
simplifies the presentation. The fluid and beam equations are coupled through the source term 
t p(u,p , h) in flHJ, which corresponds to the trace of the second component of a(u,p)ndl transported 
in the beam reference configuration. The coupling term writes: 

(4) p, h)(x, t) = —e 2 ■ a{u,p)(x,h(x,t),t){-d x h(x,t) e 1 + e 2 ), (x,t) e (0 ,L) x (0,T), 

where (ei, e 2 ) denotes the canonical basis of K 2 . The fluid and the beam are coupled also through 
the kinematic condition 


(5) u(x,h(x,t),t) = d t h(x,t)e 2 , (x,t) € (0,L) x (0,T). 

We complement our system with the following conditions on the remaining boundaries of the 
container: 

• L-periodicity w.r.t. x for the fluid and the beam; 

• no-slip boundary conditions on the bottom of the fluid container: 

(6) u(x, 0,t) = 0 . 


Note that since the question of contact is mostly a local one we assume periodic boundary con¬ 
ditions at the oulet for the fluid and replace the standard assumptions of “clamped” arterial wall 
also by periodic boundary conditions for the beam. 

An important remark on this coupled system is that the incompressibility condition together 
with boundary conditions imply: 

r L 

(7) 


I dth(x,t)dx = 0 , Vf>0. 

J 0 


Consequently, for any classical solution (u,h,p) to this system, the right-hand side of Q must 
have zero mean: 

f L 

/ (f>(u,p, h)dx = 0. 

Jo 

It can be achieved thanks to a good choice of the constant normalizing the pressure which is 
consequently uniquely defined. The pressure can be then decomposed, for instance, as 

( 8 ) p = p 0 + c, 

where, one chooses to impose 


(9) 


and where c satisfies 


>nt) 


p 0 (x,y)dxdy = 0, 


( 10 ) 


=(<) = j 


e 2 • (a(u,p 0 ))(x, h(x,t),t)(-d x h(x,t) ei + e 2 )dx. 
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This constant c is the Lagrange multiplier associated with the constraint 0. Another mathemat¬ 
ical way to compel the solution with this compatibility condition, without defining the physical 
average of the fluid pressure c, is to introduce the L 2 -projection operator on the set of L-periodic 
functions whose averages are equal to zero on (0, L), denoted M a , and rewrite 0 as 

(11) p s d t th (3d X xh T otdxxxxh r ydtxxh — M s cj)(u : p , h ). 

This is the choice made in [39] . 


2.1. Main result. In what follows, we call (BF) the “beam-fluid” system (|Tj)-([2|)-0-(|4|)-0-(j6|)- 
©-©-©■ We study herein the (BF) system, completed with initial conditions: 


( 12 ) 

h(x, 0 ) 

= h°{x), 

x £ (0, L), 

(13) 

d t h(x, 0 ) 

= h°(x), 

x £ (0, L ), 

(14) 

u(x,y, 0 ) 

= u°(x,y ), 

(x, y) £ {x £ (0, L),y £ (0, h 0 (x))} =: T° 


Our aim is to study this Cauchy problem and, specifically, to prove the existence of a unique 
global-in-time strong solution. 


We first give some notations and definitions and make precise the functional framework. For 
any given non-negative function b £ Cj {(0 , L), i.e. the set of continuous and L-periodic functions 
on R, we define 

fib := {(x, y) £ R 2 , such that x £ (0, L), y £ (0, 6 (x))} . 

With this definition F(t) = flh(t,-)- We denote by Cf°(fl 6 ) the restriction on flb of L-periodic 
functions in x indefinitely differentiable on 

O b = {(z, y) G R 2 , s.t. , y £ ( 0 , b{x))} 


Note that Ob + Lei C Ob and Ob = Ukei,fl b + Lke\. We introduce the classical spaces L\{flb) and 
H™{flb) respectively as the closures of C^°{flb) in L p (flb) or H m (flb). We define in the same way 
CJ(0,L), Lj(0,L) and H™(0,L). More generally, the subscript j) stands for the periodic version 
in the first variable of a function space. We emphasize that contrary to the usual convention, we 
consider that time is the last variable of a function. This enables to write a unified definition for 
periodic functions whether they depend on one space variable only (such as the height function 
h) or two space variables (such as the velocity-field u). Finally, we set: 


L #,o(°. L ) '■= I / e L\ ( 0 , L) s.t. J f(x)dx = 0 


and, in the same way, 

L lo(,Sl b ) := j/ £ Ljj (fl b ) s.t. J f{x)dx = 0 

Then the projection operator M s , that is applied to equation 0 leading to equation (fill) , is the 
orthogonal-projector from i 2 (0,i) onto L^ o (0,L). 

We state our main result as follows 



Theorem 1. Let us consider a > 0, /3 > 0 and 7 > 0. Assume that the initial data (h°,h°,u°) 
satisfy: 

. (h°,h°) £ Hq(0,L) x Hj(0,L), 

• u° £ H\ (J-°), 

• no-slip condition is fulfilled initially: 

(15) u°(x, 0) = 0, u°(x, h°(x)) = h°(x), Vx£(0 ,L), 

• no-contact and incompressibility compatibility conditions are fulfilled initially: 

min h°(x) > 0 , [ h°(x)dx = 0 , 

xe[o ,l] 

divu 0 = 0 onJ°. 


(16) 

(17) 


0 
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Then (BF) has a unique global-in-time strong solution. 

A precise definition of what is a “strong solution” is given in Section [3] Our proof for TheoremQ] 
follows a classical scheme: local-in-time existence and uniqueness of solutions, blow-up alternative 
and a priori estimates. Local-in-time existence and uniqueness of strong solutions has already 
been tackled in |39| for clamped boundary conditions instead of periodic boundary conditions. In 
Section [3[ we explain shortly how this result can be adapted to our functional framework with 
periodic lateral boundary conditions. This construction leads to the existence of a unique maximal 
solution for any given initial data, that blows up in finite time if and only if the quantity 

(18) C(t) := sup / + / (a\d xxx h{x,t)\ 2 +'y\d tx h(x,t)\ 2 ) dx 

xe[o,L] h ( x ’t) Jo 

fL r‘h(x,t) 

+ / A<|Vit(a;,y,t)| 2 da;d2/ 

Jo Jo 

diverges in finite time (see Corollary [5]). Note that existence of weak solutions as long as the beam 
does not touch the bottom of the fluid cavity can be obtained also by adapting Eing or m to 
our setting. 

In this paper, the main novelty is the computation, for any local-in-time strong solutions to 
(BF), of a new a priori estimate on C, defined by (TTH1) . This estimate enables us to derive 
a regularity estimate valid on any given time interval (0,T). We emphasize that, in order to 
obtain these estimates, we have to assume that a > 0 and 7 > 0. This ensures first that the 
elastic boundary remains regular, second that the beam dissipates energy. Whether or not these 
assumptions can be dropped remains a completly open question (notice that existence of weak 
solutions before contact is still valid for a = 7 = 0 and /3 > 0). 


2.2. Formal argument. Before considering the full Navier-Stokes/beam system, and in order to 
illustrate the different steps of the proof, let us first consider a reduced model for which we derive 
similar bounds (at a formal level for conciseness). This coupled system writes 

(19) Psdttb fid xx b T ol d xxxx b 7 ^xxtb — q, on (0, L), 

(20) d t b = d x [b 3 d x q\ , on (0, L ), 


where b stands for the deformation of the beam (analogue of h) and q denotes the fluid pressure. 
To complete these equations, we require that q satisfies 

(21) f q(x,t)dx = 0 Vf> 0, 

Jo 

and we impose periodic boundary conditions in x. We note that this system is related to (BF) as 
it can be obtained formally by setting h = eb and letting e go to zero. The second equation (l20l) . 
i.e. the Reynolds equation, is a reduced model for the Navier-Stokes equations (completed with 
no-slip boundary conditions) in the thin film (or lubrication) regime. We refer to [361 Section 5.B] 
for a detailed derivation of the Reynolds equation. 


Let (6, q) be a (classical) L-periodic solution to (fl9ll (1211) on (0, T) with T > 0. First, multiplying 
m by dtp and combining with (l 20 ll multiplied by q yields: 


1 d 

2 dt 


(p s \d t b\ 2 + /3\d x b\ 2 + a\d xx b\ 2 ) 


(l\dtxb\ 2 + b 3 \d x q\ 2 ) = 0 . 


We obtain that there exists a constant Co depending only on initial data for which: 


( 22 ) 


sup {p s \\d t b- L 2 (0,L)|| 2 +a||6;IL 2 (0,L)|| 2 +/3||6;IL B 1 (0,L)|| 2 ) 

te(o,T) 

+ ^ ( 7 || dtb-, H}(0,L)f + \\bid x q ; L 2 (0, L)\\ 2 ) < C 0 . 
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In what follows, Co denotes a constant depending only on the initial data but which may vary 
between lines. To obtain a lower bound on b , we multiply (1191) by — d xx b and integrate over (0, L ). 
We then integrate by parts in space. By taking into account the periodic boundary conditions 
and ( 1201 ) . we obtain: 


d_ 

dt 


Q \ d xx b \ 2 - psdtbdxxb 


+ 


(/ 3 \ d X xb \ 2 + a \ d X xxb \ 2 - p s \ d tx b \ 2 ) 


j qdxxb = j d x qd x b = J b 3 d x q d x b ) 


Finally, we deduce: 


(23) 


d t 


1 


l \ d xx b \ 2 


b 3 d x q d x 


Psdtb dxxb 


1 f L d t b 

4 p- 


'0 


{ P \ d xx b \ 2 + a \ d xxx b \ 2 ) = f p s \ d tx b \ 2 . 

Jo 


, we get that there exists a constant Co such that: 

T 


Combining with the previous bound 

(24) sup ( 7 || 6 ;^ 2 ( 0 ,L)|| 2 +||&- 1 ;C # 1 ( 0 ,L)||)+ T a\\b ; H 3 ( 0 , L )\\ 2 + p\\b ; H 2 ( 0 , L )\\ 2 < C 0 . 
tefo.T) J o 


Note that, if p s yl 0, to obtain the previous estimate one should compute an upper bound for 


pT f'L pL 

/ / p s \ dtxb \ 2 , sup / p s d t bd xx b . 

Jo Jo te(o.T) Jo 


te(o,T) Jo 

From (l22l) these terms are bounded if 7 > 0 and if 7 or a is strictly positive respectively. At this 
point, we call a real-analysis lemma which states that there exists a continuous function D m i n for 
which: 

II&- 1 ; LJ°(0,L)|| < D min (\\ b ; H 2 ( 0 , L )\\, H&" 1 ; Lj(0, L)||) 

(see Appendix IA.II for a proof). The first consequence of this inequality is that (l24l) implies that 
b remains away from zero uniformly in time. Combining this information with the dissipation 
estimate ( 1 ^ 1 ) we obtain that there exists a constant Co for which: 

sup (|r 1 (-,t);L^(0 ,L)||)+ ( T \\ q] Hi ^ L )\\ 2 < C 0 . 
t€(0,T) Jo 


Consequently the pressure is bounded in L 2 (0, T; i?j-(0, L )). Next we derive a regularity esti¬ 
mate for the deformation b by multiplying (1191) by — dtxxb . This yields after integration by parts 
in space: 


Id 

2 df 



Ps\d tx b \ 2 + a \ d xxx b \ 2 + / 3 \ d xx b \ 2 




qdt xx b . 


Thanks to the L 2 (0, T; /^(O, L))-bound on q , we reach the required estimate that enables us to 
extend solutions globally in time: 


(25) sup ( a \\ b ; H 3 ( 0 : L )\\ 2 + P \\ b - H 2 ^ L )\\ 2 + p s \\ d t b ; H ^ O , L )\\ 2 ) 
te(o,r) 

+ 7/ \\ dtb -, H 2 ( 0 , L )\\ 2 < C o . 

Jo 

This ends the formal proof of a no collision result, on the one hand, and of a global-in-time 
existence of strong solutions, on the other hand. 

In Section [2 we prove that comparable estimates hold true for the complete coupled system 
(BF). First, considering the full beam/Navier-Stokes system, the analogue of (1221) corresponds 
to the (already-known) classical decay of kinetic energy. To obtain a similar estimate to (|24l) , we 
multiply dill) by — d xx h and multiply Q by a suitable extension of — d xx h . The choice of this 
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extension is a key point of the proof (see Section 14.21) . We then obtain an identity similar to (1231) 
with additional remainder terms that we bound thanks to the energy estimate and for which a 
control of h in L 2 (0, T; H^(0, L)) D L°°(0, T; U 2 (0, L)), resp. a control of d t h in L 2 (0,T; 
is needed (hence a > 0, resp. 7 > 0). The extension of the last estimate (l25l) is more involved. 
Indeed, when dealing with the full Navier-Stokes/beam system, we also have to control the fluid 
velocity-field u in L°°(0, T ; H] (Tit))) (and not only the pressure field as for the toy model). When 
working in cylindrical domains, such an estimate is obtained by multiplying © with dtu and by 
applying elliptic estimates for the Stokes system in order to bound the convective terms. However, 
these elliptic estimates are classically proven in C ll,1 -domains or W 2 ’°°-domains mug. Here as h 
is merely L°°(0, T; H 2 (0, L)), we cannot directly apply these standard regularity properties. So, 
we need to extend the elliptic results for the Stokes system to domains which are only subgraphs 
of U 2 -functions and analyze precisely the dependency of the associated elliptic estimates with 
respect to the norms of h (see Lemma [ 6 j). This proof is an adaptation to a periodic framework of 
a lemma that can be found in [22] ■ Moreover, as the fluid domain is moving with time, instead 
of dtu, we need to consider a multiplier that takes into account this motion. The most natural 
choice is dtu + u • Vu. But, this function is not divergence-free and consequently pressure terms 
appear that cannot be handled easily. To avoid this difficulty, we mimic the method used in [9; 
in the framework of fluid/solid interactions. We introduce a divergence-free multiplier avoiding 
the introduction of the pressure in the regularity estimate. Moreorever this multiplier is chosen so 
that the associated multiplier for the structure equation is d u h. Nevertheless, a special attention 
needs to be paid since the structure motion is, once again, less regular than when considering 
fluid/solid interactions. 

The outline of this paper is as follows. In next section, we focus on the change of variables 
turning the beam/fluid system into a quasilinear system in a fixed geometry. We recall the 
construction of local-in-time strong solutions of [39| and adapt this result to our periodic boundary 
conditions framework. We end the section by a technical proposition dealing with elliptic estimates 
for the inhomogeneous Stokes system in H 2 - subgraph domains. The third and last sections 
are devoted to the proof of Theorem [T] They are divided into three subsections corresponding 
respectively to the extension of the three estimates (1221) , (l2dl) and (1251) to solutions of the coupled 
problem (BF). 


3. LOCAL-IN-TIME STRONG SOLUTIONS AND TECHNICAL LEMMAS. 

In this section, we first adapt the construction of local-in-time strong solutions of [39] to our 
periodic setting. To this end, we will apply the following change of variables: 

(26) f(x,z) = f{x,h{x)z), V (x, z) € fli. 

To measure the regularity of such a change of variable, the following technical proposition is 
required: 

Proposition 2. Let us consider h £ H 2 (0, L) satisfying min a; g[o,L] h(x) > 0. Then for any given 

m < 2 , 

• the mapping f f defined by (1261) realizes a linear homeomorphism from H™(Llh) onto 

am). 

• there exists a non decreasing function K ^ : [0,+oo[— > (0,oo) such that, if we assume 
moreover that \\h]H 2 (f),L)\\ + ||/i _1 ;Tj°(0,L)|| < Rq then there holds: 

||/;J?J"(n 1 )|| < KM)\\f-Hf(fl h )\\ , \\f-H^(n h )\\ < K^(Ro)\\f-,Hp(Sh)\\ . 

Proof. The proof is standard. For m £ {0,1} the result easily derives from the fact that h £ 
IF 1,oo (0,L) is bounded from below by a strictly positive constant. For m = 2, the key point is 
that the motion of the upper boundary is tranverse only so that we combine the regularity of h 
with the following tensorization of the space H} (lb) : 

^(Lb) = JTjfttO, L) x (0,1)) = L\ (0, L; H\ 0,1)) n H]{ 0, L; L 2 ( 0,1)). 
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The most delicate point enters the computation of \\d xx f; L 2 (f2i)||. We have: 

d xx f = d xx f + 2h'zd xy f + h'zdyf + {tiz) 2 d yy f , 

in which the worst term is h"zd v f. It is bounded in L 2 (f2i) since h" G L 2 (0, L) and 

dj e ^ l^(o,l-l 2 (o, i)) n l 2 (o,l-,l°°(o, i)). 

□ 


3.1. Construction of local-in-time solutions. As explained previously, local-in-time existence 
and uniqueness of strong solutions to (BF) are tackled in [39] with no normalizing condition for 
the pressure ((f8])-(f9]l- (fTT)ll is replaced with HU), with homogeneous Dirichlet boundary conditions 
for the fluid velocity on the part of the boundary that is not elastic and with “clamped” boundary 
conditions for the structure. Namely, instead of periodic boundary conditions, the displacement 
rj = h — 1 satisfies 


77 ( 0 , t ) = r)(L, t ) = <9^77(0, t) = d x rj{L , t) = 0 , Vi G (0, T). 

The proof of existence of solutions follows a classical method, also introduced in [2BJ [T5] when 
dealing with fluid/solid interactions. To look for solutions on a time-interval (0,T), new unkowns 
( u,p ) are first introduced applying the transformation (l26l) : 


(27) 


u(x, y,t) = u I x 


’ h(x, t)' 


p(x,y,t)=p(x 


h(x, t) ’ 


(x,y) G F{t). 


These new velocity-field and pressure {u,p) are defined in the cylindrical domain f2i x (0, T) and 
(u,p, h) is a solution to ©-©-©-©-©-©-(HU if and only if the triplet (u, p , h) is solution to a 
coupled system of quasilinear pdes that we choose not to write here for the sake of conciseness. The 
core of the existence and uniqueness result is the study of this nonlinear system. First, the author 
analyzes, via a semi-group approach, the resolution of the linear system obtained by linearizing 
around rj = 0 (or h = 1), u = 0 ,p = 0. This study is based on an accurate treatment of the added 
mass effect of the fluid on the structure through an appropriate splitting of the fluid load. Then, 
the nonlinear terms are estimated and the author proves that they remain small for a small time. 
The local-in-time existence and uniqueness of a solution to the system of nonlinear pdes is finally 
obtained by a standard fixed point argument. 


In our periodic framework, computation of nonlinearities might be reproduced without change 
while the semi-group approach might be adapted in the spirit of [35] . Consequently, for any initial 
data such that 


(28) 


h°€H$(0,L), h°€Hl{0,L), u°£H;(F°), 


and satisfying the compatibility conditions: 


(29) 

min h°(x) > 0 , [ 

h°(x)dx = 0, 


®e[o,i] J 0 


(30) 

divu° = 0, on F° , 


(31) 

u°{x,h°(x)) = h°(x)e 2 , 

x G (0, L ), 

(32) 

u°(x, 0) = 0, 

x G (0, L ), 


we obtain local-in-time existence and uniqueness of a strong solution {u,p, h ) to the Cauchy prob¬ 
lem associated with the translation of 0-©-©-®-©-@-CEU in a fixed geometry, completed 
with periodic boundary conditions. The solution verifies: 

(33) u G H\ 0, T; L 2 ^)) n C([0, T]; i/^r)) C L 2 (0, T; H 2 (n 1 )), 

(34) p G Z/ 2 (0, T; Hj[Cli)), 

(35) h G H 2 {0, T; Ljj (0, L)) fl L 2 (0, T; F7j 4 (0, L)) , /i -1 G L°°((0,T) x (0,L)) 
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We emphasize that, following the proof of [30, the pressure p is defined up to a constant for now. 
The regularity statement (l35l) ensures that the function h is Lipschitz on [0, T] x (0, L). Moreover, 
since h satisfies also 

(36) i. g W 1,oo ((0,T) x (0,1/)), 

we obtain that the domain F(t) and the non cylindrical domain defined by: 

Qt := {{x,y,s), i£(0 ,L), te(0,i), y G (0, h(s, a:))} , Vf < T, 

are both Lipschitz open subsets of R 2 and R 3 respectively. 

Going back to the moving domain by inverting the transformation (12711 . we define strong solu¬ 
tions of (BF) as follows 

Definition 3. Let the initial data ( h°,h°,u °) G H^(0,L) x H^(0, L) x satisfy (EU1) - (I3J(J1) - 

mm and let T > 0. A strong solution to (BF) on (0, T), associated with the initial data 
(h°,h°,u 0 ), is a quadruplet ( h,u,p,c) satisfying: 

• h, it, p and c have the following regularity: 

(37) h G 7L 2 (0,T;L 2 (0,L))nL 2 (0,T;iL tt 4 (0,L)), h~ x G L°°(0,T) x (0,L)), 

(38) uGHl(Q T ), V 2 uGL 2 (Q t ), cGL 2 ( 0,T), Vp € L\{Q t ) , 

• equations 0-0 are satisfied a.e. in Qt , 

• equations 0 - 0 - 02 ]) are satisfied a.e. in ( 0 ,T), 

• equations 0 - 0-0 are satisfied a.e. in (0,T) x (0,L), 

• equations (1121) - (1131) - (1141) are satisfied a.e. in (0 ,L) and F° . 

We emphasize that the pressure p in our definition is completely fixed and that u is a time-space 
function. Hence, the condition u G H^(Qt) involves both time and space derivatives of it, whereas 
Vp involves space derivatives only. The construction that we describe above, adapted from 0J], 
yields the following existence and uniqueness theorem: 

Theorem 4. Let us consider a > 0, /? > 0 and 7 > 0. Assume that the initial data (h °, h°,u°) 
belong to iJ 3 (0, L) x 47^(0, L) x and satisfy the compatibility conditions (l^l) - (15Ul) - (l3ll) - 

(13^1) . There exists Tq > 0 such that for any 0 < T < T 0 , there exists a unique strong solution to 
(BF) on (0, T). 

Proof. In the whole proof initial data (h°, h°, u°) are fixed. The only points that we want to make 
clear here are 

• the link between the regularity (I33l) - (l34l) - (l35l) of the solution (u,p,h) to the nonlinear 
system in a fixed geometry and the regularity statements of Definition [3j 

• the computation of the Lagrange multiplier c that we introduce here and that does not 
appear in [39] [38] . 

Let ( u,p , h) be the solution to III) -(I2l)-(l3l)-(l4l)-(l5l)-(l6l)- (fill) written in a fixed geometry, that one 
constructs adapting the arguments of [39 1 138 ] as explained in introduction. The deformation h(t, •) 
and its inverse l/h(t,-) are then uniformly (w.r.t. t G [0,T]) bounded in H^(0,L) and LJ°(0,L). 
So, we construct (n,p) via (1261) and apply Proposition [2] to obtain a fluid velocity and a fluid 
pressure that satisfy their contribution to (1381) . Noting that, in the aforementioned construction, 
p is defined up to a time-dependent constant and that 4 >(u,p + c,h) = 4>(u,p , h) + c, we fix p by 
requiring further that: 

(39) f <f>(u,p, h)dx = 0 , Vf e (0,T). 

Jo 

We then write p = Po + c with po satisfying 0 and with c being fixed by CQD. Due to the 
regularity of h, u and po we obtain finally that c belongs to L 2 ( 0 ,T). 
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Conversely, for any given strong solution (u,p,h,c) of (BF) in the sense of Definition (3j we 
construct (u,p) by (1371) and refer to Propostion [5] again yielding that: 

u £ nL 2 (0,T-,H$(Sh)), p £ T 2 (0,T; Hj^)). 

We apply then Bl Theorem 3.1] and deduce that u £ C([0, T); H^(Cli)) and get that, for T small 
enough, ( ii,p,h ) is the unique solution to (fTT)-0-(l^-(l4l)-(l5l)-(l6l)-(ITT1l written in a fixed geometry, 
as constructed by adapting the arguments of [39] [38] ■ □ 

Remark 1. From the regularity we just derived for u, we deduce that, for any strong solution 
(u,p,h) the mapping 1 K > fjr^ |Vrz.| 2 (t) belongs to C°([0,T]). 

Finally, we obtain that (BF) is wellposed locally in time. Following [35], it appears that we 
might choose the time Tq in Theorem [|] to be fixed by 

||/ 1 °; J ff # 3 (0, J L)|| + ||ft° ; Hl(0,L)\\ + \\h~ 1 ; Lf(0,L)\\ + ||u°; J ff # 1 (^°)|| 

only (see the computation of To at item (*), page 408). Then the following blow-up alternative 
can be classically stated: 

Corollary 5. Let a > 0, fj > 0 and 7 > 0 be given. Assume that the initial data ( h°,h°,u° ) be¬ 
long to i? 3 (0,T) x 7Jj(0, L) x and satisfy the compatibility conditions (l29l) - (l30l) - (l3Tl) - (l32l) . 

Then (BF) completed with initial conditions (I12ll - (I14|) has a unique non-extendable strong solution 
(T*,(u,p,h,c)). Furthermore, we have the following alternative: 

(■ i) either T* = +00 
(ii) either T* < 00 and 

hmsup||h(-,t);R tt 3 (0,T)|| + ||a t / l (-,t);R # 1 (0,T)|| + ||/i- 1 (-,t);Tf(0,T)|| + ||u(-,t);R tt 1 (T(t))|| =+00. 

i->T* 

The aim of Section [I] is to prove that the second alternative (ii) never holds and consequently 
that the solution is defined on any finite time inetrval (0, T). But before going any further we focus 
on the elliptic regularity properties of the inhomogenous Stokes system in a subgraph domain. 


3.2. Elliptic estimate. In this subsection we derive elliptic estimates for the inhomogeneous 
Stokes problem in a domain Qh with h £ H^( 0, L) such that ft -1 G LJ°(0, L). With this regularity, 
the domain is neither C 1,1 nor W 2 ’°° and one cannot apply standard elliptic regularity results. 
Nevertheless, we take advantage here of the fact that tth is a subgraph so that the change of 
variable transforming flh into a flat domain (namely fli) can be chosen to be smooth in the 
transverse variable (see \h below). This remark enables us to extend the classical method with 
h belonging merely to F[ 2 (0,L). Such an estimate is a key argument in the derivation of the 
regularity estimates for the solution of the non linear system (BF). 

For simplicity, we fix p = 1 in this part. Let us consider source terms and a boundary condition 

if , 9 ) e x ffjKfi*)), V G R|f (0,T). 

We aim at studying the regularity properties of T-periodic (w.r.t. x ) solutions to 

(40) — Au + Vpo = /, in fife, 

(41) div u = g , in Cl h , 

completed with boundary conditions: 

(42) u(x,h(x)) = 7)(x)e 2 , ViG (0,1), 

(43) u(x, 0) = 0, Vx G (0, T). 

Integrating div u = g over implies that the boundary velocity 77 has to satisfy 


r L r 

i)(x)dx = / 


( 44 ) 


L 


g(x,y)dxdy. 
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The left-hand side of (TCT) does not involve the deformation h, because the deformation as well 
as the boundary velocity are vertical. In what follows, we restrict to data g £ 0 (flh) and 
77 £ Ljj O (0, L) for which (Hill is clearly satisfied. 

Remark 2. Note that, for this inhomogeneous Stokes problem, with Dirichlet boundary conditions, 
the pressure po is defined up to a constant. Consequently, we enforce uniqueness of the pressure 
by imposing: 

/ Po(x,y)dxdy = 0. 

■Nth 

The main result of this section writes 


Lemma 6. For any h £ H^(0,L) such that h 1 £ source terms and boundary condition 

C f,g) e L^su 0 x (^(fi*) nLjf i0 (fi h )), f, e ff*(o,L)ni£ 0 (o,L), 

there exists a unique solution (u,po) £ x nLj 0 (fi/,)) to the Stokes system (|40)l - 

(HU-(H2J)-(I43D. Moreover, there exists a non-decreasing function K s : [0,oo) —> (0,oo) such that, 
if we assume \\h; H^(0, L)\\ + ||h _1 ; L|°(0, L)|| < Ro then, this solution satisfies: 

(45) \\u-,H^n h )\\ + Ibo; ^(14)11 < K s (Ro ) (||/;Lj(n fc )|| + + lb; h}(0, L)\\) . 

The remainder of this section is devoted to the proof of Lemma [6] This proof is an adaptation, 
to our periodic framework, of the computations on the Stokes problem that can be found in [22] . 
which, itself, uses ideas of [IS]. Compared to [22] , we also carefully analyze the dependance of the 
constant K on h in inequality (1451) . To obtain the expected dependency, we assume throughout 
this section that: 

|b;Lr # 2 (0,L)|| + ||/ i - 1 ;Lf°(0,L)|| < R 0 

and show that K depends only on Rq. 


First step: Rewriting of the Stokes system in a given geometry. As in [22] we compute regularity 
estimates for solutions to (l40l) (l43l) by studying the Stokes system transported in a geometry which 
does not depend on the deformation h. Namely, we derive regularity estimates on (u,p) defined by 
(EH. Indeed, thanks to Proposition [21 we remark that (u,p) £ H^(VLh) x (f lh) is a solution to 
(l40ll (l43l) if and only if (u,p) £ H^(VLi) x H^(Q 1 ) is a solution to the following Stokes-like system 

(46) - div[(AhV)u] + (BhV)po = /, in Hi , 

(47) di y(B^u) = g, in Hi, 
completed with boundary conditions: 


(48) 

(49) 


u(x, 1 ) = f](x)e 2 , \/x£(0,L), 

u(x, 0) = 0, Va;G(0,L), 


where ( Ah,Bh ) and (/,<?) are explicit. Indeed, by introducing the mapping Xh{x,z) = (x,h(x)z), 
for (x,z) £ Hi, we obtain: 


(50) B h := cof V X h = (J f , A h := i(cof S7 Xh ) T cof V Xh 



-h'z \ 

I 1 (.h'z ) 2 I > 

h ~ r h ) 


and the transported source terms: 

/ := hf, g := hg. 

We note that / £ L 2 (H 1 ), g £ Lj^ 0 (fii) D iJ^fii). Thus, thanks to Proposition [2] to prove 
Lemma[ 6 ]it is sufficient to derive similar estimates but on the transported unknowns (u,p) solution 
of ([4(1 [49]) . Namely, we obtain that there exists a unique (u,po) £ R 2 (H 1 ) x (iJ^H) nL 2 0 (Hi)) 
solution of C 6 l) ( 1121 ) that satisfies 

(51) Ib ;j ff # 2 (fii)|| + ||po;H„ 1 (n 1 )|| < j a(||/; j L 2 (H 1 )|| + ||s;£^i)|| + \\rpH 2 (0,L)||) , 
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where the constant K depends only on R 0 . Since the matrices Bh and Ah are in 7/ ((0, L); H s ( 0,1)), 
for any s > 0, we have that Ah and Bh belong to a multiplier space of H 1 ( f2i). We refer the 
reader to [221 Lemma 6] for more details. In particular, we obtain that for any v £ H 2 (VLi), there 
holds div[(A/jV)i;] £ L 2 (Lli) and, for any q £ there holds (B;,V)q £ Thanks to 

Piola identity, we also have: 

(52) di v(Bjv)=BZ :Vv, 

so that, for any v £ H 2 (Ll i), there holds div(B^v) £ Consequently the assumptions on 

the deformation h are compatible with the expected regularity on (u,p). 


Remark 3. Let us mention that we get estimates for a pressure po such that 

/ po(x, z)dxdz = 0 . 

Through the change of variables (12611 , this implies that the pressure q, defined by q( x, y ) = 
Po{x,y/h(x)) and on which we deduce an estimate, verifies the following constraint 

f q(x,y) . 

Jn h K x ) 


-dxdy = 0. 


So, the pressure we compute with this method does not match the one mentioned in Lemma [21 
Nevertheless, the effective pressure po mentioned in LemmaUH reads Phq where Ph stands for the 
L 2 -orthogonal projector on Li 0 (Clh): 

Po = P h q :=q- 


|fb, 


q(x,y)dxdy 


satisfying ||po; H^(flh)\\ < ||g; Rj 1 (J7/ l )||. Hence, we prove (1451) with po replaced by q. 


Let us now study precisely the existence and uniqueness of (u,p) in Hl{ f2i) x L 2 Q (Lli) and 
derive elliptic estimates in H 2 {TLi) x H^( f2i). 

3 

Second step: Lifting of the Dirichlet boundary conditions. As f] £ H£(0,L) there exists u^ £ 
H 2 (Vl i) such that uy | _ i = rje 2 , u^, _ 0 = 0 and 

(53) IK;tfj?(fii)|| <C\\t,-,H*( 0 ,L)\\. 

with a constant C depending only on the fixed geometry. We set u = u — u^. This new velocity 
satisfies: 


(54) - div[(A,,V)u] + (B h V)p 0 = f, in Oi, 

(55) div(Rju) = g, in fli, 
with 

f ■= f + div[(A,,.V)u^], g := g - Bj : Vuf, , 
completed with boundary conditions: 


(56) u(x, 1) = 0, Vx £ (0,L) , 

(57) u(x, 0) = 0, Vx £ (0, L). 

As underlined previously, thanks to the regularity of Ah and Bh, the new sources terms ( f,g ) 
belong to L 2 (fli) x HffTli) and satisfy the following estimates: 


(58) II/; L \(nOH + ||g; J ff # 1 (fIi)|| < tf(||/;£?(fii)|| + \\~g; ^(^i)!! + \\y-,Hf (0,L)\\), 

where K depends only on Rq. Moreover the average of g on fli is still equal to zero, since 



Recalling that Ur, satisfies we obtain now that the proof of Lemma [6] reduces to the study of 
the case rj = 0 ( i.e. solving system d5Tl) - (l55l) - (l5Bll - (lF7ll '). 
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Third step: H 1 x L 2 estimates. We first define as the dual space of the subset of (O-j) 

of functions with zero trace on (0, L) x {0} and (0, L ) x {1}. The aim of this step is to prove that, 
for any (f,g) £ x 0 (fii), there exists a unique (u,po) £ -^(fli) x 0 (fli) solution of 

([5i | - (f57| and satisfying 

(59) + ||po;£?(n 1 )|| < K (||/;H # - 1 (Q 1 )|| + ||ff;£j(ni)||) , 

where if depends only on Rq. Since the arguments are quite standard (see, for instance, [25], [221 
in similar contexts), we only sketch the main points of the proof. First, we notice that: 

• Ah £ L°°(Qi) and there exists two non negative constants a\ and ai controlled by above 
and from below by a function of Rq for which 

ail < Ah(x, z) < c^I) V (x, z) £ flj , 

in the sense of symmetric matrices; 

• Bh is invertible and Bf ^ 1 belongs to ifj(( 0 , L)\H s ( 0 , 1 )), for any s > 0 , with norms 
dominated by a function of Rq only. 

With the second point at-hand, we build a lifting operator for the divergence. Namely, for any 
X £ L 2 0 (f2i), there exists a function w £ Hj(fli), with w \ z=1 = w\ z=0 = 0, such that 

(60) div (B^w) = x, ||ni:F # 1 (n 1 )||<if||x;^(^ 1 )||, 

where K depends only on Rq. Indeed, as x has zero average on fli, there exists v £ Hj( fli), with 
t)i i = ui 0 = 0 , such that 

div (v) = x, ||i> : tf||(fti)|| < C||x;i#(^i)||- 

See, for instance, na Lemma III.3.1]. We set then w = B h T u. As B h T is a multiplier of H 1 with 
norm bounded by a function of Rq we obtain (lC>Ul) . 

Then, to solve (lMl) - (fo7| we first lift the divergence source term g by applying the previous 
construction. We then solve the Stokes-like system (I54l) - (l55l) by reproducing the classical argu¬ 
ments for the Stokes system. As Ah satisfies the first point, we first construct a weak solution 
u £ By(fli) depending continuously on ( f,g ). Then, as B ^ 1 satisfies the second point, we obtain 
also the pressure p 0 £ L 2 0 (fli) which completes (l59l) . 

Fourth step: proof of Lemma\fh H 2 /H 1 -regularity. To complete the proof of Lemma[ 6 l it remains 
to obtain an estimate on the second order derivatives of u and the first order derivatives of pq. We 
obtain that 

(61) \\u,H 2 (Sl i)|| + ||po;-ff B 1 (fI 1 )|| < K(\\f-L 2 {Sh)\\ + ||g;ir # 1 (fii)||) , 

where K depends only on Rq. We follow the method introduced in [¥S] and already applied in 
[22] in our subgraph framework. Thanks to a classical regularization argument, we assume in 
what follows that h £ Cfi°(Q,L). In this case, classical elliptic estimates ensure that (u,po) £ 
ff 2 (fli) x i/j (Oi). Nervertheless the standard elliptic estimates involve norms of the deformation 
in W 2 , 00 ( 0, L) (see for instance 0). Consequenlty, we aim to show that the constant only involves 

Ro- 

First we obtain estimate on u x := d x u and p x := d x f>Q. For this purpose, we differentiate the 
equations (l54l) . (l55l) satisfied by (u,po) w.r.t. x. We obtain that ( u x ,p x ) £ fli) x L 2 0 (fli) is 
the solution of 

—div[(A^V)u x ] + ( B h X7)p x = f x , on Oi, 

di v (B^u x ) = d x g — d x B^ : Vu, on fli, 

where f x = d x f + div[(9 x A^V)?l] — (d x Bh'V)p o , completed with periodic boundary conditions on 
lateral boundaries of fli and homogeneous boundary conditions on y = 1 and y = 0 (we recall 
that we consider the case 77 = 0 ). 

We note that 


d x Bj : Vu = di v(d x B^u), 
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which implies 



: Vu = 0 . 


Consequently, taking into account that g is L-periodic w.r.t. x, we obtain that g x = d x g — d x B^ : 
Vu has a zero average on Oi. Due to the regularity of (u,po), the right hand side (f x , g x ) belongs 
to 7Jj _1 (f2 1 ) x Ly.(Sli), but we need sharp estimate to show our main result. As we stated previously 
(■ d x Ah,d x Bh ) € L 2 ((0, L ),; H s (0, 1)), for arbitrary s > 0, (with norms bounded by a function of 
Ro) and Hj(( 0, L ) x (0,1)) C Lj°((0, L);L 2 ( 0,1)). Hence, in the spirit of [22] Lemma 6], we obtain 


(62) ||div ((a c A, l V)iz);H # - 1 (n 1 )|| < ||(0 x 4 h V)«;Lj(n 1 )|| < K\\u~, ff # 1 (fi 1 )|| 1/2 ||u*; F # 1 (0 1 )|| 1 /2 ) 
and 


(63) || d x B h : V«;Lj(fi 1 )|| < iL # 1 (D 1 )|| 1 / 2 ||ii x ;iL # 1 (D 1 )|| 1 / 2 , 

where if depends on i?o- Next we have to estimate (d x BhS7)po in (Qi). Thanks to the Piola 
identity and the fact that Bh is the cofactor matrix of the gradient of Xh , we obtain, for any 
w € TJ^(Hi) such that tO| i=0 = w\ z=1 = 0 

/ (d x B h V)p 0 w = - Pod x Bj : Vw. 

Consequently, as in the computations of the latter bounds, we obtain: 

(64) ||(^ J B h V)p 0 ;i? # - 1 (fii)|| < ^||p 0 ;L 2 (H 1 )|| 1 / 2 ||p x ;L 2 (H 1 )|| 1 / 2 . 

We can now apply the result obtained at the previous step to (u x ,p x ). Combining with (l62l) . (l63l) . 
HMD, this leads to 

||A c ;iL # 1 (L! 1 )|| + Wp^LK^W < I< (||/;Lj|(fi 1 )|| + \\g; H'(n^W + ||u; 7L # 1 (D 1 )|| 1 / 2 ||u x ; ^ 1 (H 1 )|| 1 / 2 

+ ||po;Lf(Di)|| 1/2 ||p x ;L^(Hi)|| 1/2 |) . 

and finally: 

(65) \\u x ; i?j 1 (fl 1 )|| + \\p x ; L 2 (H 1 )|| < K (||/; L 2 (H 1 )|| + \\g- HjM) . 


To obtain a similar estimate on the full second order gradient of u (resp. on the full gradient 
of po), we have to bound d zz u (resp. d z po). To this end, we note that, differentiating (l55l) w.r.t z 
and applying (li)5l) . we have 

II - zh'd„ui +d zz u 2 -,L 2 (n 1 )\\ < iL(||/;L 2 (D 1 )|| + ||s;Ljf(fii)||) . 

While, combining the first equation of (l54l) with the second equation of (1541) multiplied by zh', in 
order to eliminate the pressure, leads to 

\\ztid zz Ui+d zz u^q{Sh)\\ KKiWf-L^W + Wg-H^n^W). 

By simple algebraic combinations, since 1 + z 2 (h') 2 > 0, we obtain 

IIV 2 ^; Lf (fii)|| < K (||/; ^(D 1 )|| + ||£; Ljf(fii)||) . 

A similar inequality holds for || d z p\ L 2 (Hi)||, using once again the first equation of (1M1) . Combining 
these inequalities, we finally obtain the desired bound 

||«;iL 2 (B 1 )|| + ||po; J ff # 1 (^i)|| < K (||/;Lj(f2 1 )|| + ||ff;^(D 1 )||) . 

This ends the proof of Lemma [Gj 

For the study of the whole coupled system, we need an estimate on the surface load applied by 
the fluid on the structure. If we compute 4>{u,p , h) through the change of variable (l27l) (be. w.r.t. 
u,p and h,) and we use, for instance, the multiplier Lemma EH Proposition B.l] or Proposition 
[2] we can also obtain the following corollary, stated without proof: 
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Corollary 7. Let h £ H^(0,L) such that h 1 £ L^°(0,L) be given there exists a non decreasing 
function K b : [0,oo) —> (0,oo) such that, if \\h', H^(0, L)\\ + ||h _1 ; Lj°(0, L)|| < Rq the following 
propositions hold true. 

Given source terms ( f,g) £ L^(£lh) x D L 2 0 (£lh)) and a boundary condition i) £ 

3 

H£(0,L) (~l Lj o (0,L) satisfying (l44l) . the unique pair (u,po) £ H^(£th) x (H^(Qh) D L^ 0 (£lh)) 
solution to m - nTD - (l^ - C3l) satisfies: 

(66) \\(f)(u,po,h)-,H^(0,L)\\ < K b (R 0 ) (||/;Lj(n h )|| + ||<?; ^(ft/OII + ||r); ^ § (0, L)||) . 

The constant c defined by c = £ J { ^ f>{u,po, h)dx satisfies 

(67) |c| < K b (R 0 ) (\\f;Ll(n h )\\ + \\g; H # 1 (^)ll + \\t h} (0, L)||) . 


4. Proof of Theorem []] 

Let (h°,h°, u°) £ R$(0, L) x Hj( 0, L) x be given and satisfy the compatibility conditions 

(l29l) - (l32l) . We consider (u,p, h, c), the associated non-extendable strong solution to (BF) completed 
with the initial conditions (fl2ll - (fl4l) (in the sense of Definition [3]). This solution is defined on some 
time-interval [0,T*), where T* > 0. We compute estimates satisfied by this solution on [0,T] for 
arbitrary T <T*. 

The proof of Theorem[l]is divided into three parts, each of them corresponding to the derivation 
of one estimate similar to (l22l) . (l24l) and (l25l) respectively. First, we recall the energy estimate 
satisfied by the solution, then we prove a distance estimate which ensures that the beam does not 
touch the bottom of the fluid cavity on the time interval (0, T) and finally we derive a regularity 
estimate which garantees that the strong solution can be extended on any given time interval, 
leading to our global-in-time existence theorem. 


4.1. Energy estimate. We first recall the classical estimate associated with the dissipative equa¬ 
tions that we consider. We introduce £ c and R respectively the total energy of the coupled system 
and the dissipated energy: 


£c (t) := \ 


(p s \d t h(x,t )\ 2 + a\d xx h(x,t )\ 2 + j3\d x h(x, t)\ 2 ) dz + 

U{t) := 7 f \d tx h(x,t)\ 2 dx + p I \X7u(x, y, t)\ 2 dxdy . 

Jo Jr(t) 


'o 
■L 




p f \u(x,y,t)\ 2 dxdy 


>Ht) 

We have then that: 

Proposition 8. The following energy balance holds true 

( 68 ) 


£ c (t)+ f H(s)ds = £ c (0 ), Vi € [0,T]. 

Jo 

Proof. Multiplying first m by u and integrating by parts leads to 

n (dfU + u ■ V«) ■ u = / / divcr(u,p)-u 

7 (s) Jo Jf(s) 

= f [ a(u,p)n-u- 2 p f [ |-D(m )| 5 

Jo JdJ r {s) Jo Jf(s) 

pt pL pt 

(69) 


/ 0 JdT(s) 
r t pL 


[ [ <t>{u,p,h)d t h- 2p [ j \D(u)\ 2 . 
Jo Jo Jo Jj-(s) 


0 JT(s) 

In the last equality we have used the coupling conditions at the interface between the fluid and the 
structure. Moreover thanks to the only vertical motion of the beam together with the divergence 
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free constrain (see [3; Lemma 6] in the 3D case), there holds: 


(70) 


2 [ \D(u)\ 2 = [ |V 

Jt(s) 


Furthermore, since the boundaries of J-(t) move with the velocity-field u , we have 

S=t 

r x r i r 

(71) 


[ (dtu + u ■ V«) ■ u = — 

/ l«l 2 

h(s) 2 

Jr{s) 


s=0 


Consequenlty 


\u°\ 2 - 


/ / </>(u,p,h)d t h. 

J 0 Jo 


(72) \f \u\ 2 + pj f |Vu| 2 = i/ 

z j j-{t) Jo Jr(s) z Jt 0 

If we now multiply the beam equation ([3]) by dth we obtain, after time and space integration by 
parts, 


(73) 


1 

2 



{ps\dth\ 2 + a\d xx h\~ + P\d x h\ 2 ) 


S=t 


s—0 



4 >(u,p, h)d t h . 


o Jo 


By summing (1721) and (1731) . we obtain the expected result. 


□ 


4.2. Distance estimate. The aim of this section is to prove the following proposition 


Proposition 9. There exists a constant Co depending only on initial data for which: 

r-T 

sup ( 7 ||h(t,-);77 # 2 (0,L)|| 2 +||h- 1 (t,-);C # 1 (0,L)||)+a / \\h(t, •); i7 # 3 (0, i)|| 2 df < C 0 (l + T). 
te(o,T) v 7 Jo 


The remainder of this paragraph is devoted to the proof of this result. Let us consider 0 < t < T. 
In all what follows Co denotes a constant depending only on the initial data but which may change 
between lines. Let w = = (—d y tp, d x ijj) where: 

y, t) = d x h(x, t)x o > V (*> *) e Qt , 

with: 

Xo(z) = z 2 (3-2z), \/z 6(0,1). 

Combining the regularity of h (which implies that h £ C([0, T]-H^(0, L)) D iJ 1 (0, T; L7 2 (0, L)) 

thanks to [3DJ Theorem 3.1] see (1371) 1 with xo 6 C' oo ([0,l]) we obtain that w £ L7 1 (Qt) and 

V 2 w £ L 2 (Qt). This regularity is enough to justify all computations below as ( u,p,h,c ) is a 
strong solution. Moreover, w is divergence free by construction and, since Xo(l) = 0 and X'o(l) = 
Xo(0) = Xo(0) = 0, there holds: 

w(x,h(t,x),t) = d xx h(t,x)e 2 , for all x £ (0, L) and t £ (0, T), 

w(x,0,t) = 0, for all x £ (0, L) and t £ (0, T). 


Consequently, we multiply m by w and © by d xx h and integrate on Q t for arbitrary t < T. 
We get after integration by parts (note that the terms involving the fluid/beam interactions cancel 
out since the structure test function is the trace of the fluid test function on the interface): 


f pf(d t u + u-\7u)-w-2fj J [ D(u) : D(w) + f f (/3\d xx h\ 2 + a\d xxx h\ 2 ) 
Jqa. Jqa. Jo Jo 


(74) 


+ 


( 2 \ ®xxh\ PstJthdxxh' 


t r L 


J s =o 


Ps | &tx h | • 


0 j 0 
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We first show that this identity leads to an estimate that is comparable to (1731) up to remainder 
terms we shall bound afterwards. The term that will enable us to bound /i -1 is 2y Jg D(u ) : D(w). 
To deal with this term, we introduce a well chosen pressure 

(75) q(x,y,t):=q s (x,t) + d xy ip{x,y,t), q s {x,t):=- d yyy ip(s,y,t)ds, V (x, y, t) € Qt ■ 

Jo 

d x h(x t') 

An easy computation gives d yyy il>(s,y,t) = — 12yA—-T—, so that 9s satisfies 


(76) q s (x,t) = 12f 

J 0 


d x h(s,t) 


ds = 6 


( h(x,t )) a 
1 1 


0 (MM)) 3 [\h( 0 ,t )\ 2 \h(x,t )\ 2 _ 


V(x,f) g (0, L) x (0, T). 


In particuliar, q s does not depend on y. Furthermore Vg G T 2 (Qt)- Applying again the fact that 
w is divergence-free, we obtain: 


(77) 


! [ D{u) : D(w) = [ (2 D(w) - ql 2 ) ■ D(u ) 

JQt JQt 

(2D(w) — 9 / 2)71 • u — / (Aw — Vg) ■ u . 
) J Q.± 


0 JdF{s) 


Note, that all the terms make sense thanks to the regularity of (w,q). In this last identity, by 
definition (1751) of q, we have: 

/ (Aw-Vg)-w = / d xxx ip u 2 - 2d yxx i\) m 

JQt JQt 

n (n 1 d xx 'i/ju 2 - 2n 2 d xx ipu 1 ) da - d xx ip(d x u 2 - 2d y u 1 ) 

'JF(s) JQ t 

= - / d xx ip(x,h(x,s),s)d t h(x,s)d x h(x,s)dxds - / d xx i/j(d x u 2 - 2d y ui). 

Jo Jo J Qt 

Similarly, the other term of (1771) writes 
[ f (2D(w) — ql 2 )n ■ u 

JQ JdJ-(s) 

= JJ d t h(x,t)((d yy ip(x,h(x, s),s) - d xx ip(x, h(x, s),s))d x h(x, s) 

+ d yx il>{x, h(x, s ), s) - q s (x , s)^dxds . 

Differentiating the identity d y i/j(x, h(x, s), s) = 0 (that holds true since Xq( 1) = 0), with respect 
to x, yields 

d yx ij)(x, h(x, s),s) + d x h(x, s)d yy ijj(x, h(x, s), s) = 0. 

Consequently, we simplify: 


n (2D(w)-qh)n-u=— / d t h(x,t)d xx ip(x,h(x,s),s))d x h(x,s)dxds 

>JF(s) Jo Jo 

/ / dth(x,s)q s {x,s)dxds. 

Jo Jo 

Combining the computations of both terms in (1771) . we obtain finally: 

(78) 2 f D(u) : D(w) = - f [ d t h(x,s)q s (x,s)dxds+ f d xx ij){d x u 2 - 2d y u{) ■ 

J Qt Jo Jo J Q t 
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At this point, we replace q s by its explicit value (see (1761) 1 and by remembering that the average 
of dth is zero, we obtain 


(79) f f dth(x, s)q s (x, s)da;ds = — 6 f f 

Jo Jo Jo Jo 


t r L 


d t h(x, s ) 
/o Jo \h(x, s)| 2 

Consequently, from (l78l) . (1791) . the equality (l74l) reduces to: 


deeds = 


6 


h(x, s ) 


da: 


s=0 


(80) 


J (jj\d xx h\ 2 - p s d t hd xx h+ +J J (/3\d xx h\ 2 + a\d xxx h\ 2 ) 

/ J s=0 

= / / Ps\d tx h\ 2 +n / d xx i /)( d x U2 - 2 dyUi ) + / Pf{d t u + u • Vu) • w 

Jo Jo JQt JQ t 

We recognize in the left-hand side of this equality the quantities that we want to estimate as in 
Compared to (l23l) . we have two additional terms 


Ti = p d xx tl){d x u 2 - 2d y u {), T 2 = pf(d t u +u ■ Vu) ■ w. 

JQt J Qt 

To bound these terms, we need precise estimates on the stream-function ip that are gathered in 
Appendix IA.2I 

First, Ti is bounded by applying Proposition fl5l and energy estimate 


'Qt 


(Jxx'lp(JJx'U j 2 ^dyll l) 


< C\\d xx iP; L 2 (Qt)\\\\Vu-, L 2 (Q t )\\ 


< 


Co [||/i; ij°(0, L )|| || d xxx h; L 2 ( 0, L)|| 5 


+|| dxxh; L$( 0, L)p\\d xxx h- L]{ 0, L)||i 

with a constant Cq depending only on the initial data. Once again, using the energy estimate 
we obtain 




< Co 


sup ||/ l ;LJ°(0,L)||5) ( f \\d xxx h;L 2 (0,L)\\' 
.te(o,T) / Vo 


+ 1 sup^ \\d xx h] Ljj (0, A)|| 4 ] ( j \\d xxx h; L^O, L)|p 


(81) 


< C'q(1 + T)+£i / \\3 xxx k;L 2 i{ 0,L W 


for arbitrary small £\ > 0, that will be choosen later on. We note that Cq depends on ey a priori 
but the value of this parameter will be fixed to a universal constant afterwards. This remark is 
also valid when other e’s that are introduced. 

Concerning T 2 , taking into account the convection of the fluid domain by the fluid velocity, we 
have, 


-| S=t 

(82) T 2 = (dtu + u ■ Vu) ■ w = j u(-,s) ■ w(-,s) — j (dtiv + u-'Vw) 

JQt _Jf(s) J s=q Jq ± 


u. 


We bound the first term on the right-hand side of (|82|) using a Cauchy-Schwarz inequality. Apply¬ 
ing the energy estimate (l68l) and the estimates (1122I) - (I128I) on the gradient of the stream-function, 
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one gets 


r 

s—t 



L 

s =0 


< C 0 + \\u(-,t);L 2 (Tm\\m-:t):L 2 (Tm 


(83) 


< C 0 I 1 


< Cq + £2 


dx 


h(x, t) 


dx 


h(x,t) ’ 


for arbitrary small £2 > 0 . 

For the second term of the right hand side of (15^1) . we first integrate by parts in space. Since 
i/>(x, h(x, s), s) = d x h(x,s) and d y ip{x, h{x, s), s) = 0 , we have dti/j(x,h(x,s),s) = dt x h(x : s) so 
that: 

pt pL 

u • dtw = — 

iQt JO JO JQ t 

rt r L 


n dth(x,s)d x h(x,s)dtt/j(x,h(x,s),s)dxds + / {d y u\ — d x U 2 )dttp 

JQt 

n d t h(x, s)d x h(x , s)d tx h(x , s)dxds + / ( d y ui - d x U 2 )d t 4 >. 

JQt 


Applying the energy estimate (1681) and using the following ID embedding inequality 

\\d x h;L^(0,L)\\<C\\d xx h;L 2 (0,L)\\, 
we can estimate the boundary term 


n d t h(x, s)d x h(x , s)d tx h(x, s)dxds 


< 


II d t h- L\{ 0, -Zb)|| || d x h- Lf (0, -Zb)|| || d tx h- L\{ 0, L)|| 

J 0 

< c f \\d xx h-,L 2 (o,L)\\\\d tx h-,L 2 (o,L)f 


(84) < Co ■ 

Taking into account (I124D and the energy estimate (1551) . we now bound the second term by 

< C 0 \\dtip; L 2 (Q t )\\ 


(■ d y ui - d x u 2 )d t 4> 


(\\dt x h;L 2 t (0,L)\\ 2 + \\d xxx h;Ll(0,L)\\) 


< Co 

(85) < Co(1+T)+s 3 [ ||a rax / l ;-Lj(0,L)|| 2 , 

J 0 

where £3 > 0 will be chosen later on. 

Finally for the last term in the right hand side of (|82|) . we have, after space integration by 
parts: 

r t r L 

( 86 ) 


/ u-Vu>-u= / / \dth{x, s)\ 2 d xx h(x, s)dxds — / u-S7u 

JQ t Jo Jo JQ t 


Concerning the boundary integral in (1861) . we apply (1681) to show the following estimate 


n \d t h(x, s)\ 2 d xx h(x, s)dxds 


( 87 ) 


< C (\\d xx h(-, s); Lj ( 0 , 1 /)|| \\dth(-, s); ( 0 , T)|| 2 ) ds 

Jo 

< sup \\d xx h; T«(0, A)|| [ \\d tx h; Li (0, L )|| 2 < C 0 ■ 

te(o,T) Jo 
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Moreover for the volume integral in the right hand side of (1861) . we have 
u ■ Vm • w 


'Qt 


< 


t fL I / rh(x,s) \ 2 / ph(x,s) \ 2 

/ hi 2 / |W| 2 sup \w{x,y,s)\ | dzds. 

0 -'O \ WO / Wo / j/e( 0 ,/i(a,«)) 


From (11211) we know that the following pointwise estimate holds true 


\w(x,y, s)| < C ( \d xx h(x,s)\ + 


\d x h(x,s)\ | d x h(x,s)[‘ 


Thus we define 


t nL [ f nh(x,t ) 


and 


h = 


Ii = 


h = 


h(x, s ) h(x, s ) 


rh(x,s) \ - \ 

J |Vm | 2 j \d xx h(x, s)| I dzds, 


\/{x,y) € J'(s). 


0 ^0 \ Wo 


* rLI/ r"\uA 

o Jo \ \ Jo J \J o ) h(x,s) 


da;ds 


t fL / ( ph(x,s) 


0 JO \ \J 0 




dxds. 


We now take care of each quantity. Applying the ID embedding inequality 

\\d xx h-,L^(0,L)\\ < C\\d xxx h-,L 2 (0,L)\\ 
and the energy estimate (1S51) . we obtain: 


h < C 


< C 0 


o \ 

•t 

'0 \ \JF(s) 


) (^ () l Vu ! 2 ) \\d xx h;L^(0,L)\\jds 

j[[J l Vu n \\dxxxh\L}{G,L)\\ J ds 


< Coll Vu; Li(Q T )\\ 


\\d x xxh> (0, L)\\ ds 


< ——b £4 [ \\d xxx h' : L 2 (0, L)|| 2 ds, 

£4 Jo 


for arbitrary small £4 > 0. We now take care of 1 2 . We note then that u vanishes on y = 0 so that 
a Poincare inequality yields: 

/ ph(x,s) \ 5 ( Mx,s) \ 5 

(/ M 2 J < Ch(x, s) ( |Vu| 2 J . 

Using this bound to estimate I 2 we get 

h < C f \\\/u-q{T{s))\\ 2 \\d x h-L^(Q,L)\\As<Co 
Jo 

as \\d x h; LJ°(0, L)\\ < C\\d xx h; T|(0, L)\\ which remains uniformly bounded in time (see (l68l) l. 
With similar arguments, we also prove I 3 < Co- This yields finally 


u ■ V ii • w 


'Qt 


< Cq + £4 J \\d xxx h m , (0, U)|| 2 ds 
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and thus, taking into account (1571) 


( 88 ) 


u-Vwu 


< Cq + £4 / \\d xxx h; £jt( 0 , C)|| 2 ds . 


Finally, T 2 can be bounded, thanks to (1531) . (1571) . (1551) . (1551) . as 

(89) |T 2 | < C 0 (l + T) + £ 5 T || d xxx h- L \{ 0, L )|| 2 + £ 2 t -^-r 

Jo Jo ri\x,t) 

where £5 and £ 2 are to be chosen small enough. 


Combining (1591) and (15T1) to bound the right-hand side of (l80l) and taking into account (1651) to 
bound the remaining terms on the right-hand side depending on h that might be concerned, we 
get, for any t <T, 



;\d xx h(-,t)\' 


_ 6 P_\ 


+ 


[ [ (/3\d xx h\ 2 + a\d xxx h\ 2 ) 
Jo Jo 


< Co(l + T) + PsWdM; i); L\i 0, L) IIII d xx h(-, t)-qi 0, L) || 

+Ps j J M+eU hi^t) + j ^ xxx ^^i 

< Co(l + T) + e ^ ^ + / ||^;^(0,L)H 2 d S 

for some arbitrary small £ > 0. We conclude the proof of the distance estimate of Proposition [9] 
by choosing e small enough. 

Remark 4. 

• The derived distance estimate relies strongly on the fact that the beam motion is only trans¬ 
verse and that we control the curvature of the elastic boundary. Indeed, we need to have 
bounds on the deformation h in L°°(0, T; R 2 (0, L)) and in L 2 (0, T; H^(0, L)) to control 
remainder terms. Both norms are bounded because a > 0, even if the first one is controlled 
via the energy bound while the second one is controlled simultaneously with h~ l . 



• To prove our distance estimate we need to control dtp in L 2 (0, T; i?j(0, Lf). This is the 
reason why we assume 7 > 0. One may wonder whether the fluid dissipation would be 
sufficient. A priori, from the L 2 (0, T; H^IFft)) bound of the fluid velocity we only get a 

control on dth in L 2 (0, T; H^ 2 (0, L)), which is not enough. 

4.3. Regularity estimate. Combining Proposition [9] with Proposition H71 we obtain that h^ 1 e 
L°°(0, T; Lj(0, L)). More precisely, we have that there exists a non-decreasing function R : 
[0,T*) —y [0, 00) bounded on all bounded subintervals of [0, T*) such that: 

(90) ||h- 1 (-,t);L oo (0,L)|| + ||h(-,t);R 2 (0,L)|| < R t Vte[0,T*). 

In particular, the beam never touches the bottom of the fluid cavity on bounded time intervals. 
This lower bound on h ensures that the elliptic regularity result derived in Section 13.21 applies 
and enables us to pursue further in order to prove global existence of strong solutions by choosing 
appropriate test functions. We obtain the following quantitative estimate: 

Proposition 10. There exists a function C° eg : [0, T*) —> [0, 00) bounded on all bounded subin¬ 
tervals of[0,T*) such that: 

(91) ||u(qt);H„ 1 (^(t))|| 2 + ||ftM^);^ 1 (o,i)ll 2 + IIM^);^^o ) L)|| 2 <c^(t) 

The remainder of this subsection is devoted to the proof of this proposition. We fix T < T* 
and construct the C° eg for t £ [0,T]. We split the proof into three steps. In the first one, we aim 
at multiplying the fluid equation and the structure equation by dtu and dtth respectively. Yet, 
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d t u is not an appropriate multipier of the fluid equation since it does not take into account the 
motion of the fluid domain. A natural choice is then the total derivative dtu + u ■ Vit, but this 
function is not divergence free. So we introduce a modified divergence free test function following 
ideas of [9;. This step requires to bound the fluid velocity in H 2 (P(t)) and consequently, thanks 

to the elliptic estimates derived in Section lT2l the structure velocity in H^ 2 (0,L). As we do 
not get estimates on this quantity in our first step, we need a second step to obtain a regularity 
estimate on the deformation of the beam. This second estimate depends itself on the regularity of 
the applied fluid-force and thus on high-order norms of the fluid-velocity and the pressure. The 
final step consists in a well chosen combination of the two previous estimates in order to obtain 
the expected result. 

Remark 5. 

• As the function R : [0,oo) (0, oo) satisfies (IflUl) at-hand, when we apply below Proposi¬ 

tion [H Lemma[6\ or Corollary the associated respective constants K e ,K s or K b define 
non-decreasing functions of time bounded on bounded subintervals o/[0,T*). 

• In the computations below, we denote by Cq a constant depending on inital data and by 
Co : [0,T*) —> (0, oo) a function that is bounded on all bounded subintervals of[0,T*). The 
value of this constant/function may vary between lines. 

Step 1 : L°°(0,T; H 1 ) regularity estimate of the fluid and structure velocities. As explained in the 
previous paragraph, we introduce a vector-field A(-,f) which coincides with u on dJ-(t) for all 
t £ (0,T). and which satisfies 


• A is divergence free, 


( 

. vA e L 2 /Qt) and v 2 A e L 2 (Qt) 

with 


(92) 

||A(.,i);Ljf(.F(i))|| 

< 

C o (t)\\d t h(;tf,L 2 (0,L)\\, 

(93) 

||vA(.,t);L 2 (^(t))|| 

< 

Co (t) || d t h(-,t) ;H}(0, L )||, 

(94) 

||V 2 A(.,t);L 2 (J-(t))|| 

< 

C o mdth(;t)-H 2 (0,L)\\, 


for a.e. t £ (0, T), 




• A = u and c^A = 0 on y = h{x, t ), 

• A = 0 on y = 0. 

The construction of A is given in Appendix [5] With the notations of this appendix, we have 
Co(t) = K l (R t ) that is indeed a function which is bounded on all bounded subintervals of [0,T*). 

Next we define v as 

v := d t u + A ■ Vu — u ■ VA. 

Given t < T, it is a suitable multiplier for (HD on Q t as it belongs to Lj(Q t ). Indeed, by classical 
Sobolev embedding, we have 

IIA • Vu; T 2 (Qt)|| < / t ||A;i 4 (^))||||V U ;L 4 (^( S ))|| 

Jo 

< C 0 {t) / t ||A;^ 1 (-A(s))|||k;R 2 (J-( S ))|| 

Jo 

< C 0 (t) sup \\d t h(;s);H;(0,L)\\ (|| W; L 2 (Q t )|| + || V 2 u; L 2 (Q t )ll) . 

se(o,t) 

We note here that we have used the continuous embedding L 4 (J r (s)) H 1 (J 7 (s)). A priori, 

the constant associated with this embedding depends on the domain JF(s) and thus on the de¬ 
formation of the beam. However, going back in a fixed domain and interpolating the results of 
Proposition [2] (or extrapolating an equivalent version for the L 4 -space that we skip for concise¬ 
ness), we might prove that this constant is uniformly bounded locally in time, as || h{-, s); L)\\ 
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and ||h _1 (-, s); L|°(0, L)\\ are bounded locally. In the same way, in what follows, we may use in¬ 
terpolation inequalities in JF{s) for which the constant will be bounded by a locally bounded 
function of ||h(-, s); 17 2 (0, L )|| and ||h _1 (-, s); LJ°(0, L)\\ and thus by a locally bounded function of 
R s . Similarly, we obtain u-VAeI 2 (Q ( ). 

So, we multiply m by v on Q t and obtain the following identity 

(95) / pf(dtu + u ■ Vu) • (dtu + A • Vu — u ■ vA) = / diver • (dtu + A ■ Vu — u ■ VA). 

Jr(t) Jr(t) 

Formally, the trace of v on y = h is equal to dtth. So we expect that this identity has to be 
combined by the structure equation multiplied by dtth. Again, dtth is a suitable multiplier for (0 
as it belongs to Ljj((0,L) x (0, t)) and we obtain the following identity: 

pL pL pL pL pL 

(96) / p s \d tt h\ 2 — / cj>(u,p,h)d t th =-a d xxxx hd tt h + 3 I d xx hd tt h + 7/ d xxt hd tt h 

Jo Jo Jo Jo Jo 

To compute the right-hand sides of (1551) and (ITOl) . we need more regularity than the one satisfied 
by the strong solution under consideration. Yet the following lemma holds true: 

Lemma 11. For any triplet ( w,q,b ) satisfying the regularity assumptions of Definition^ namely 

h G H 2 { 0, T; L 2 { 0, L)) n L 2 (0, T; ff # 4 (0, L)), 


w G , V 2 w G L 2 (Qt) i g 6 L^(Qt) , Vq G L 2 (Qt) , 

and such that w(x,h(x,t),t) = dt.b(x,t) G L^ o (0, L) and divw = 0, the following identities are 
satisfied: 


t r-L 


(97) / divcr(u), q) ■ (d t w + A • Vw — w ■ VA) = — 


'Q t 

%[ \Vw\ 2 + w / |Vui| 2 (0) — 2p f D{w) : ([vA] T Vu; + vA[Vw] T — D(w ■ VA)') , 

2 JT(t) 2 Jr( 0) JQ t v 2 


(j)(w,q,h)d tt b 


0 jo 


and 


pZ pL/ pZ p L-j 

dxxxxbd tt b — fi / / d xx bd tt b-'y / / d xxt bd u b = 

Jo Jo Jo Jo 

pL 

^ Q \d tx b(x,t)\ 2 - pd t b(x,t)d xx b(x,t) - ad tx b(x,t)d xxx b(x,t)j dx 

n L pt pL 

\d tx b\ 2 a / \d txx b\ 2 

Jo Jo 

pL 

- ^ (J^\d tx b(x,0)\ 2 - pd t b(x,0)d xx b(x,0) - ad tx b(x,0)d xxx b(x,0)^j dx. 

The proof of this lemma relies on regularization arguments and is postponed to Appendix [B] 
Thus, we apply (l97l) for w = u and b = h. Then identity (l95l) becomes 


t r L 


(99) / pf(dtu + u ■ Vu) ■ (dtu + A • Vu — u • VA) = — 


4 >{u,p, h)d tt h 


0 .70 


h 


IVul- 




£ f |Vu 0 | 2 -2p [ D(u) : ( [VA] T Vu + VA[Vu] T - D(u ■ VA)) , 
2 Jr{ 0) JQt V ' 
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For the left-hand side of (1M1) . denoted LHS, we have: 


LHS = 


/ pf\d t u + u ■ Vu\ 2 — / Pfiptu + u ■ Vu) • (u • Vu) 

JQt J Qt 


+ pf (dtu + u ■ Vu) • (A ■ Vu — u • VA), 


i j 


St 


p/|d t u + u ■ Vu\ 2 — — / p/|u • Vu| 2 + / p/ (dtu + u ■ Vu) • (A • Vu — u ■ VA), 


'Qt 


which yields 

(100) ^ f |Vu| 2 + i f p f \d t u + u- Vu\ 2 + [ [ <j>(u,p,h)d tt h 

z JF(t) z JQ t Jo Jo 

<— f pf\u-S7u\ 2 — j pf(dtu + u • Vu)(AVu — u ■ VA) 

2 dc t JQt 

-2 p D{u) : ( [VA] T Vu + VA[Vu] T — D{u ■ VA)] + w / |Vu 0 | 2 . 

dQt v / 2 7^(o) 

We split the right-hand side of this inequality into six integrals denoted /ithat we bound 
independently. Applying interpolation inequalities for estimating the L 4 -norm, we have 


/i := 


Pf\u • Vu | 2 


< C? 0 (i) / ||n;L 4 (J-( S ))|| 2 ||Vn;L 4 (J-( S ))|| 2 d S 

Jo 

< Co(t) / t || M ;L 2 (J-( S ))||||Vn;L 2 (J-( S ))|| 2 ||u; J ff 2 (J-( S ))||d S . 

Next we use the elliptic estimates derived in section [3T2l to bound ||u; i^ 2 (J 7 (s))||. 

(101) ||u; i? # 2 (J-(s))|| < X s (i? t )(l|5 t u + U • Vu; Ljf(.F(a))|| + ||d t h; JT* (0, L)||). 

Following Remark 0 the function of time t K s (R t ) is (non-decreasing) and bounded on all 
bounded subintervals of [0,T*). Thus, there holds: 


h < 


Co(t) f ||«; L] (.F( a ))||||Vu; L 2 (T(s)) || 2 (|| d t u + u ■ Vu; L 2 (T(s))\\ + \\d t h; h} (0, L)||) da. 

J 0 

Recalling the energy estimate proven in Proposition [8] we obtain, for arbitrary e > 0, 

(102) h < C 0 (t) f ||Vu; A 2 (J r (s))|| 4 ds 
Jo 

+ e [ (II d t u + u • Vu; Ll(F(s ))|| 2 + || d txx h; Lj( 0, L)|| 2 ) ds . 

Jo 

Here and in what follows, we skip the dependance of Co(t ) on e as this parameter will be fixed 
universally later on. The second term reads: 

h'■= / Pf(d t u + u ■ Vu) ■ (A • Vu), 


that we estimate as follows 

rt 


< 


\\dtu + u ■ Vu■ L 2 (F(s ))\\II A; L 4 (J-( S ))|| || Vu; Lf(.F(a))||da 


< C 0 (t) \\d t u + u-Vu-,q(T(s))\\\\Vu-,Ll(R( s ))\\i... 

Jo 

• ■ • l|u; ^ 2 (-? r (s))|| ^ || A; I/ 2 (J r (s))|| 3 1 | A; ^jj 1 (^ 7 (s))|| ^ds . 
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Thanks to estimates (1^21) . (l93l) satisfied by A and using Proposition [8] which enables us to bound 
||u; L^(F(s))\\ and \\dth; A 2 (0, L)\\, and the elliptic estimate (11011) . we get, for arbitrary e > 0 to 
be fixed later on: 

(103) \h\<C 0 (t) C (|| Vu; Lj(J r (s))|| 4 + \\d tx h; A 2 (0, A)|| 4 ) ds 
Jo 

: [ (||a tM + U -Vu;Lj(^( S ))|| 2 + ||5 tra / l ;Lj(0,L)|| 2 )d S . 

Jo 


The third term reads: 


We estimate it as follows 


+ e 


h ■— / Pf(d t u + u ■ Vm) • (u ■ VA). 
J Qt 


I3 I < 


II d t u + u • Vu; A 2 (J'(s))||||u; A 4 (.F(s))|||| VA; Ljf(.F(s))||ds 


< C 0 (t) [\\d t u + u-Vu-,Ll(T(s))\\\\u-,q{ns))\\^ ■■■ 

Jo 

... || Vu; Ljf(.F(s))||* || A; J ff # 1 (J'( S ))|| * ||A; tf # 2 (.F(s))|| ^ds, 

Applying estimates (l93ll - (f94|l satisfied by A, together with interpolation inequalities and Proposi¬ 
tion [Si we obtain 

|/ 3 | < C 0 (t) f || d t u + u- Vu; L\(JF(s)) || || Vu; L\(T(s)) || * || dt x h; ij( 0, L)|| *1| d txx h; L\{ 0, L)|| *ds.. 
Jo 

Hence, for arbitrary e > 0, there holds: 

(104) |/ 3 |<C7o(t) /' i (||Vu;^(J'( S ))|| 4 + ||a tx / l ;ij(0,A)|| 4 )d s 

Jo 

+ e [ (||d t u + u- Vu;^(J-(s))|| 2 + ||a txx / 1 ;^(0,L)|| 2 )d S . 

Jo 

We proceed with 

h ■= [ D{u) : ([VA] T Vu + vA[Vu] T ), 

jQt 

that we bound as follows (C is a constant depending only on p): 

\h\ < C / t ||Vu; J L^( S ))|| 2 ||vA;L 2 (J-( S ))||d S 
Jo 

< Co(t) /* l|Vu;i^(J'(s))||||u;iJj 2 (J : '(s))||||A;7J # 1 (J r (s))||ds. 

Jo 

Similarly as above, using (f93l) and (|101D . we obtain, for arbitrary 5 > 0: 

(105) \I 4 \<C 0 (t) f (\\dtxh; Ajj (0, A)|| 4 + || Vu; A 2 (.F(s))|| 4 ) ds 


+ £ 


f (II d t u + u ■ Vu; A^(.F(s))|| 2 + || d txx h- Aj(0, A)|| 2 ) ds . 
Jo 


Finally, the fifth term is defined by 

h ■= / D(u) : D{u • VA). 


Expanding D(u ■ VA), we obtain 
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The first term on the right-hand side is bounded as I 4 (see (11051) ). As for the second term of the 
right-hand side, we have 

/ |Vu||u||V 2 A| < C 0 ® f ||A; A? # 2 (^-( S ))||||z i ; J L2(^-( s ))||i||v^; J L2(_ 7r(s ))|||| xt ; jff 2(^( s) )||i ds 


< 


J 0 

Consequently, we obtain: 

/ |Vm||u||V 2 A| < C 0 (t) f ||V«; L 2 (T®)\\ 4 
jQt Jo 

+ e [ (\\d t u + u-\/u® 2 (R(s))\\ 2 + \\d txx h® 2 (0,L)\\ 2 ) . 
Jo 

Finally, for arbitrary small e > 0, we have: 

(106) \I 5 \<C 0 ® C (\\dt x h] L^(0, A)|| 4 + ||Vu; L 2 (J r (s))|| 4 ) ds 


+ £ 


f (||$u + u ■ V«; L 2 (T(s)) II 2 + \\dtxxh; L 2 { 0, L)|| 2 ) ds . 

Jo 


Introducing (ll02l) - (fT051) - (ll0dl) - (ll05D - (11061) into (11001) . we obtain that, for arbitrary e > 0, there 
holds: 


(107) 77 / |Vu | 2 + j [ pf\d t u + u ■ Vu | 2 + f [ <j)(u,p,h)dtth 

z Jr(t) 4 JQ t Jo Jo 

< Co® f (||Vu;L 2 (J r (s ))|| 4 + ||5 te / i ;L 2 (0,L)|| 4 ) 

Jo 

+ £ [ {\\9tu + u-'Vu;L 2 (J : (s))\\ 2 + \\dtxxh',L 2 (0,L)\\ 2 ) + ^- f |VmoT 

Jo 1 Jf( 0) 


We now take care of the elastic part: we apply (|98l) of Lemma fill for b = h and rewrite |M|) as 


(108) J J p s \dtth\ 2 + J (JJ\d tx h{x,t)\ 2 - fid t h(x,t)d xx h(x,t) - ad tx h(x,t)d xxx h(x,t)^ da; 

pt pL pt pL pt pL 

-P / \dt x h\ 2 — a / \d txx h\ 2 = / / (j>(u,p, h)d tt h + C 0 

Jo Jo Jo Jo Jo Jo 

We add (11071) and (11081) . restricted to e < min(p//8,1) and bound all possible terms by energy 
estimate (Proposition [8|. This yields: 


(109) ^ [ |Vu| 2 + i [ pf\dtu + u- Vu| 2 + / f p s \d tt h\ 2 + Rem® 

z Jr(t) ° JQt Jo Jo 

<Co® f (|| Vu; Lp(J r (s))|| 4 + \\dt x h; L 2 (0, L)|| 4 ) + (a + 1) f [ \d txx h\ 2 + C 0 . 
Jo Jo Jo 


where, thanks to Proposition [8] 


pL 

Rem® = J (Jj\d tx h{x,t)\ 2 -/3d t h(x,t)d xx h(x,t) - ad tx h{x,t)d xxx h(x,t)) dx, Vt€[0,T]. 


satisfies: 


Rem® > -Co - ^ (\\d xxx h(-,t); L 2 (0, L)\\ 2 + \\d tx h(-, i); L 2 (0, L)|| 2 ) . 
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Finally, we get: 

(110) |V «| 2 - ^ / (\d xxx h(x,t )\ 2 + \d tx h(x,t)\ 2 ) dx 

z Jr(t) z Jo 

+ ^ f Pf\d t u + u- S7u\ 2 + [ I p s \dtth\ 2 
° JQt Jo Jo 

<C 0 (t) C (|| Vu; L 2 (J r (s))|| 4 + \\dt x h; £ 2 (0, T)|| 4 ) + (a + 1) [ [ \d txx h \ 2 + C 0 . 
Jo Jo Jo 

To go further, we thus need to obtain an estimate for 


f [ L \d txx h\ 2 

Jo Jo 


{\d xxx h(x,t )\ 2 + \d tx h(x,t)\ 2 ) dx 


The next step is thus an independant regularity estimate for the beam equation. 

Step 2: Regularity estimate for structure equation. First, — dt xx h € L 2 ((0, L) x (0, T)) so it is a 
suitable multiplier for 0. We obtain after integration: 


pt nix pt plx pt plx pt pL 

/ / P s&tth*Jtxxh “b (3 I I 3 xx hdt XX h Oi I j 3 xxxx hdtxxh T y / / 

Jo Jo Jo Jo Jo Jo Jo Jo 

n 4>(u,p, h)dtxxh. 


We can further integrate by parts 

r L 


(111) i/ [p s \dt x h(x, t )\ 2 + a\d xxx h(x, t )\ 2 + /3\d xx h(x,t)\ 2 ] dx + 7 / f \d txx h \ 2 = 
z Jo Jo Jo 

-[ [ 4>(u,p, h)d txx hdx+ i [ p s \d x h 0 (x )\ 2 + a\d xxx h 0 (x )\ 2 + /3\d xx h 0 (x)\ 

Jo Jo z Jo 1 


Note that all terms make sense due to the regularity of the strong solution h. We control the first 
term of the right-hand side with the help of Corollary [7J recalling that constant K b (R t ) defines 
a function of time that remains bounded on bounded subintervals of [0, T*). Consequently, we 
obtain: 


</>(«, P, h)d txx h 


Since pAu — Vp = d t u + u ■ V« we have then 


< U{u,p,h)-Hf{Q,L)\\\\d txx h-H^ 5 (0,i)|| 

< Co(t) (\\pAu — Vp; I/ 2 (J r (t))|| + \\dth; (0, T)||) \\d t h; H?(0,L)\\. 


<j>{u,p, h)d txx h 


< C 0 (t) (llfltu + u • Vit; L 2 {F(t))\\ + II &th\ HJ ( 0 , L)||) \\d t h- Hf ( 0 , L )||. 


Moreover, we apply the following interpolation inequality 

\\dth; H£ (0, L)|| < C\\d tx h-, L 2 (0, L )\\* \\d txx h-, L 2 (0, L )\\* . 

this yields: 

r L 


(112) 


</>(u,p, h)d txx h 


< C 0 (t)\\d tx h ; Lj(0, L )|| 2 + e (|| d t u + u ■ V«; L 2 (T(t ))|| 2 + \\d txx h; L]{ 0, L)|| 2 ) . 
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Introducing (1112D into (I111D . we obtain that, for arbitrary e > 0 sufficiently small, there holds: 


-i /* nl n 

(113) - [p s \d tx h(x 7 t )\ 2 + a\d xxx h(x,t)\ 2 ] dx + j- / / \d txx h \ 2 

z Jo z Jo Jo 

<C 0 (t) [ \\d tx h; Lj( 0 , L )\\ 2 + e f \\d t u + u ■ Vu; Tj((.F ( S ))|| 2 + C 0 
Jo Jo 

Step 3 : Conclusion. We multiply (|113|) by a sufficently large constant, typically: 


max 


a 4 4 ( a + X ) 


,P 8 7 

and add (11101) . we obtain, choosing e sufficiently small: 


(114) % [ \Vu\ 2 +[ \d t u + u-Vu \ 2 

1 Jr(t) JQ t 

[p s \d tx h(x,t )\ 2 + a\d xxx h(x,t)\ 2 ] dx + ^ J J |<9 tra ft| 2 ^ 

<C 0 (t) [ (l + || Vit; L 2 (J r (s ))|| 2 + \\dt x h- £ 2 ( 0 , £)|| 2 ) 2 + Co , 
Jo 




where <5 > 0. Setting 


£ - = 
^reg ■— ^ 


J |Vu| 2 + 5 [\J [Ps\d tx h{x,t)\ 2 + a\d xxx h(x,t)\ 2 ] dxj. 


we have then for all t £ [0,T]: 

IK', t)-H\T{t))W 2 + ||d*ft(-, t); £^(0, L) || 2 + || ft; ff 3 (0, i) || < Co + C£reg(t). 
with a constant C depending only on p, p S7 a and S , and 

£reg(t) < C 0 {t) f (l + || Vu; L 2 (T(s)) || 2 + \\d tx h' 7 L^(0, L)|| 2 ) (f« fl (s) + l)ds + C 0 
Jo 

Proposition [H] implies then that: 

f (1 + ||V ? i;L 2 (J-(f ))|| 2 + ||^ft;L 2 (0,L)|| 2 ) < C 0 (t) 

Jo 

We thus complete the proof of Proposition [TO] by applying a Gronwall lemma and remembering 
that t i-a C'o(i) is bounded on any bounded interval. 

Appendix A. Technical details for the no-contact result 

In this appendix we collect technical lemmas that are used throughout Section [I] 

A.l. Estimating positivity of scalar functions. We prove first functional inequalities which 
enable to bound from below a positive function. 

Proposition 12. Let a £ (0,1/2) \ {1/4}. Given a non-negative function p £ H 2 (0,L), there 
holds: 

Jo ^ CWdxxP] -b 2 (0, L) || 2 1| 77 ; ( 0 , Z/)|| 2 *- 1_2a) . 

Proof. Without restriction, we assume that 77 is smooth and we set 

b(x) = [? y(x)] 1_ “ Vx£(0 ,L). 

A straightforward integration by parts together with the L-periodicity, yields 


|c/z&(x)| 4 dx + 3 / b(x)d xx b(x)\d x b(x)\ 2 dx = 0 . 
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Replacing b by its value, we obtain: 

(l-a) 3 (l-4a) / dx + 3(1 - a ) 3 ^ |t?(: x)\ ( - 1 ~ 2a) d xx rj(x) 


/o 


\v(x) | 


\dxV(x)\' 

\v(x)\ 2a 


■dx = 0 . 


When a ^ 1/4 this yields 

r L \dxV(x )\ 4 

J 0 


rL 


\v(x)\ 


d x<C 


4a — 


|?y(a ;)| (1 2 a) d xx r](x) 


\dxV(x)\‘ 

\ r l( x )\ 2a 


dx. 


We conclude by applying Cauchy-Schwarz inequality, which yields, since a < 1/2 



\dxV(x )\ 4 

\rj(x)\ 4a 


dx 


< 

< 


C / i | ? 7 (x)| 2 ( 1 - 2 “)|9 ra 7 ? (x)| 2 dx 
Jo 

C||r ? ;Lf°(0,L)|| 2 ( 1 - 2 “)||5 xx 7 ? ;L 2 (0,L)|| 2 . 


□ 


Proposition 13. Given a non-negative function rj £ iL 3 (0,L), the following pointwise estimates 
hold true: 

(115) \d x p(x )\ 2 < C\\d xx r); Ljj (0, L)\\i\\d xxx ip, L 2 (0, L)|| ^ |? 7 (x)|, Vx £ (0 ,L). 

(116) \d x rj(x )\ 2 <C\\d xxx r)-L^(0,L)\\\r](x)\, Vx£(0 ,L). 

Proof. Let us denote 2 £ H^(d,L) a non-negative function. 

(117) \d x z{x)\ < \\\d x z \ 2 + zd xx z-L^°(0,L)\\^ , Vx€ (0,1). 

Indeed, as z is L-periodic, d x z reaches its maximal and minimal values. Then, if d x z is maximal 
at xo £ ( 0 , L ), then d xx z(x o) = 0 and consequently 

\d x z{x 0 )\ < \/\d x z{xo )\ 2 + zd xx z(x 0 ) < \\\d x z \ 2 + L^°(0,L)\\i . 

We get a similar inequality where d x z is minimal and obtain (I117D . 

We apply now the previous estimate with 

z{x) = \fp{x ), Vx £ (0, L). 

Of course z is a non-negative L-periodic function which belongs to H^(0,L) so it satisfies (11171) . 
Replacing 2 by y/rj leads to 

\dxV(x)\ < CyJ\\d xx r]-L^(0,L)\\ r]{x ), Vx£(0,L). 

Thanks to the continuous embedding (0, L) D L 2 o (0, L) in L^°(0,L), we obtain (I116D and the 
interpolation inequality \\d xx rj\ Lj°(0, L)|| < \\d xx r]; L 2 (0, L)\\ 1 / 2 \\d xxx ry, L 2 (0, L )|| 1/2 implies that 
(I115D is satisfied. 

□ 


Proposition 14. There exists a continuous function D min : (0,oo) x R + —> (0,oo) such that, 
given a non-negative function rj £ H^(0,L) there holds 

(118) W 7 ]- 1 j ^tT(0, ^)ll — 77mm (||i? _1 ; Lj(0, L)\\, || 77, 77 2 (0, L)||) . 

Proof. Without further restriction, we assume p to be smooth. We note that || 77 ” 1 ; L“(0,L)|| is 
achieved for some xq £ [0,L]. Then, as d x rj(xo) = 0 we have: 

r](x) = rj{xo) + / (s - x 0 )d xx r](s)ds , 

J Xo 

At this point, we consider two cases. First, if \\d xx rj; H 2 (0, L)\\ = 0, we obtain 

(119) ||r ? - 1 ;Lj(0,L)|| = J L||7r 1 ;Lf(0,L)||. 
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Second, if \\d xx r]; H^(0, L)\\ > 0, we have the bound: 


r){x) < vM + -/=\\ir,H$\\\x £c 0 | 2 , Va; e (x 0 - L,x 0 + L). 

For simplicity, let denote r] m := rj(xo) and M = ||ry; U 2 (0, L) ||/3. We have then 


rXo+L i 

\\V > 

J Xn 


ds 


> 


1 


'x 0 ^( s ) J 0 r]m + M\s\i ~ Mirit 

where we applied a straightforward change of variable to introduce 

rLa dz 


M 


<K*) ■■= [ 

JO 


/ o (i + M*) 

We remark that : R + —> R is an increasing continuous function such that 

{ La for a « 1, 

/' dz 

/ -;r for a » 1. 

J R 1 + |z|5 

Applying furthermore that \\rj~ 1 ; Lj(0, L)\\ < L||ry _1 ; L“(0, L)|| = Lr ] m , we obtain 

'M||r ? - 1 ;Lj(0,L)|r 


Finally, if || 77 ; -Hj 2 (0, L)|| > 0, we have : 
( 120 ) 




||7T 1 ;L5“(0, j L)|| < ^^(O, J L)||ir,- 1 ;Lj(0,L)|| 


||7 ? ;i/ # 2 (0,L)||||r ? - 1 ;Lj(0,L)|| 


-1 -3 


The above expansion of 4>{a) for small values of a yields that, when || 77 ; H^(0, L)\\ « 1 with 
|| 77 —1 ; Lj(0, L)\\ bounded, there holds 


|| ?7 ;^(0,L)|| 2 ||7 ? - 1 ;L}(0,L)|r 

Hence, we set 


||7 ? ;g 2 (0,£)||||r;- 1 ;^(0,L)|| 


OL C? 

Dmin(a, /3) = — , if /3 = 0, D min (a , / 3 ) = —— 

-L y 


f3a 


2 \ A-3 


-||»?- 1 ;L}( 0 ) L)||. 


else. 


Because of the previous arguments, this is a continuous function on (0, 00 ) x [0, 00 ) which satisfies 

(fim □ 

A.2. Estimates on the stream-function %). In this appendix, we gather technical estimates 
regarding the stream-function 

V>(z, V , t) = d x h(x, t)x 0 ([ ) V e ® T ’ 

where h £ fJ 2 (0, T; (0, L)) fl L 2 (0,T; H^(0, L)) is given and remains strictly positive on (0, T). 
We recall that 

Xo(~) = z 2 (3 — 2z) Vze( 0,1), 

that J-(t) stands for the fluid domain at time t and Qt is the time-space fluid domain. With these 
notations, we prove: 


Proposition 15. There exists a universal constant C < 00 for which: 
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for all t G (0, T) 


(121) \V^(x,y,t)\<C 


\d xx h{x,t)\ 


\d x h(x,t)\ \d x h(x,t)\ 2 


h{x,t) h(x,t) 


(122) < C\\h-,L^(0,L)\\i\\d xx h;Ll(0,L)\\i 


\/{x,y) G T(t) 


da; 


h(x, t) 


(123) \\d x ^{;t)-L\F(t))\\ < C\\h-,Lf>(0,L)\\?\\d xx h-,L 2 (0,L)\\ ; 

• on the whole time-space domain 

r T 


(124) 


L 2 (Q t )\\ < C 


(125) 


\\d xx ip;L 2 (Q T )\\ < C 


{\\h-L^(0,L)\\\\d xt h-L 2 (0,L)\\ 2 

+ \\dt.h; -Lj (0, L)\\ 2 \\d xxx h; £ 2 (0, Z/)||^ 
£ (||^; ^J°(0, L)\\\\d xxx h; Lj (0, L)|| 2 

+ \\d xx h-L 2 (0 1 L)\\i\\d xxx h;L 2 (0,L)p 


Proof. We start by the time-dependant estimate. We compute the partial derivative of the stream- 
function with respect to the space variables: 

V \ | d x h{x,t)\ 2 y , ( y 


h(x,t) J \h(x,t)\ 2 Xo \h(x,t )J ’ \/( xv t)GOx 
d x h{x,t) , ( y ^ V{x,y,t)eti T . 


d x tp{x,y,t) = d xx h(x,t)x o 

d y if(x,y,t) = ' ' xb i , / I ) 

h(x,t) \h(x,t) J 

Thus (1121D is satisfied. Furthermore, we have, for all t G (0,T), recalling that xo £ C°°([0,1]): 

,l /■/*(*,*) |a xA(x> t)|2 


\\dyipi L 2 (T(t ))\\ 2 = [ [ 

Jo Jo 


o Jo 
L r 1 


Xo 


\h(x,t)\ 2 

\d x h(x,t)\ 2 , 2 


2/ 


Jo Jo \h(x,t)\ 

Then applying Proposition 1121 with a = 3/8, we have: 


h(x, t) 

IXo(*0r<tydz. 


dydx 


||a^;L 2 (J-(t))|| 2 < C 


f 


\d x h(x,t)[ 


■dx 


dx 


/i(x, i) 


o \h(x,t)\ 

<C Wh-Lf^L)^ 

I/O | 

<C ||/i;LJ o (0,L)||^||a xx / l ;L 2 (0,L)|| 
With similar arguments, we easily obtain 


/' L \d x h(x, t) 4 

2 

i- 

1 _ 

io |/i(x, t)| 5 


-1 

if 

o 

_ 1 


dx 

h(x,t ) 


\M\L 2 (m)f < C II h- 0, L)\\\\d xx h-L 2 ( 0, L) II 2 + £ ^f dx 
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Applying again Proposition [l^] with a = 3/8, the following estimate holds true in a similar way 
as above 

\\d x ip\ L 2 (T(t ))|| 2 < C\\h- 0, L)\\\\d xx h- L 2 { 0, L )|| 2 . 

This ends the proof of (I122I) - (U23I) . 

We now take care of the two other estimates. A simple calculation gives 

y 


dti/j(x , y , t) = d xt h(x, t)x o 
Consequently, we have 


h{x , t) 


d x h(x,t)d t h(x,t)y , ( y 
~Xo 


\h(x,t)\ 2 


, V (a;, y, t) £ Qt- 


||^;£ 2 (Qt )|| 2 <C \\h-,L£>(0,L)\\\\d xt h-,Ll(0,L)\\ r ' 


h(x, t) 

L \d x h{x,t)\ 2 \d t h(x,t)\ 2 
h(x, t ) 


da; 


With Proposition [13] we bound the ratio \d x h(x, t)\ 2 /h(t, x) in the right-hand side, and we obtain 

\\dtip;L 2 (Q T )\\ 2 < C f T (||/i; ij°(0, A)||||9 a:t /i; L 2 (0, L)\\ 2 + \\d t h\ L 2 {f), L)\\ 2 \\d xxx h; L 2 (0, L)||) . 
Jo 

Finally, we estimate d xx ip which is equal to 


d xx if(x, y, t) = d xxx h(x,t) Xo( ,/ ,J ~ 

n[x, t) 


y s 3yd xx h(x,t)d x h(x,t) , y (d x h(x,t)) 3 y 2 


Xo(T7T7T)+ 


(■ h{x,t)) 2 

Consequently, the following estimate holds true 

^ f \h-Lf{Q,L)\\\\d xxx h-L 2 (Q,L)\\ 2 

\d xx h(x,t)\ 2 \d x h(x,t)\ 2 |<9 x /i(a;,t)| e 


Xo( , / .J , 


h(x,t) ( h(x,t)) 4 ' h(x,t) 


\\d xx ip-,L 2 (Q T )\\ 2 < C [ 

Jo 


h(x,t ) 

By estimates (11151) and (11161) of Proposition [13] we obtain 
\\d xx if;L 2 (Q T )|| 2 < C [ T \\\h;L^(0,L)\\\\d xxx h;L 2 (0,L)\^ 


\h(x,t)\ 3 


■da; 


■|| d xx h\ L 2 (0, T)|| 2 ||9a:xa:^j (0, L)|| + ^d xx h\ Tj(0, L )|| 2 \\d xxx h; L 2 ( 0, L)|| : 


Finally (11251) is satisfied. 


□ 


Appendix B. Technical lemmas for the regularity estimates 

In this second Appendix, we collect technical results that are applied in the proof of Proposition 
M We first construct the lifting velocity-field A and prove then two identities that are central in 
the proof of our regularity estimates. 

B.l. Construction of A. We construct a time-frozen operator Uh which satisfies the equivalent 
properties to the one we require for A. We shall then set A = Uh[dth], The construction of Uh is 
the content of the following proposition: 

Proposition 16. Let h £ H 2 (0,L) such that h~ l £ Lj°(0,L), there exists a continuous linear 
mapping U h : H^(0,L) (~l L 2 q (Q,L) -a s.t. : 

• for all r) £ Hj( 0, L) D L 2 0 (0, L) we have 

U h [rj\(x, 0 ) = 0 , U h [ri\(x,h(x)) = r)(x)e 2 , Vx£(0,L), 
div C/ft. [ 77 ] = 0 on Llh 

• for all rj £ H 2 (Q,L) D L 2 o (0,L) there holds Uh[r)\ £ H^(VLh) . 

Furthermore we construct Uh such that: 

• d 2 U h [fj\(x,h(x)) =0. 
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• there exists a constant K l depending increasingly on \\h; H^(0, L)\\ + \\h Lj°(0, L)|| for 

which 

(126) \\U h [%L](Sl h )\\ < K l \\r)-, I/ 2 (0, L)||, 

(127) \\U h [r,Y,Hl(n h )\\ < K^rj-H 1 (0,L)||, 

(128) \\U h [fi\-,Hf(Sl h )\\ < K l ||»?; H 2 (0,L)||, Vi) £ H 2 ( 0 , L) n L£ 0 (0, L ). 

Proof. We consider h £ H?( 0, L) which does not vanish on (0, L) and assume that 

||/ l ;i7 # 2 (0,L)|| + ||/ 1 - 1 ;Lf(0,L)||<i? o . 

We construct Uh in order that (I126D - (U27D - (I128D holds with a constant K l depending on R 0 only. 
Fix A = l/(2i? 0 ) so that min^gjo.L] h(x) > 2A. Assume that 77 £ H^{ 0, L). We define Uh as 

[ rje 2 , in Cl h \ (0, L) x (0, A), 

(129) U h [r,\ = { 

[ Ux[fj\, in (0 ,L) x (0, A). 

Note that the restriction of Uh to £lh\( 0, L) x (0, A) obviously satisfies continuity estimates similar 
to (1126l) - (1127l) - (1128l) . With this step, we reduce the construction to a rectangular box. In order to 
change the boundary on which we have to match the data, we define U\[rj\ by: 

U\[v](x,y) = (^ q ^ U[rj\ > V( X ’V) e (0 ,L) x (0,A). 

Let us define now U[rf\. We expand g £ H^(0,L) in Fourier series as 

v{x) = Y Vn exp (in^-). 

n£Z 


By assumption, we have 770 = 0. Let k £ N. Then, choosing Q a polynomial function such that: 

(130) Q(0) = 1 Q'(0) = Q" (0) = 0 

(131) Q w ( 1) = 0 \/l<k + 1. 
and, for all n £ Z: 

P„( 0 ) := <5(min(|n|2,1)) Vze[0,oo), 
we obtain that P n £ C fc+1 ([0,1]). Consequently, we define 

U[r)\ = V^, with ^(a;, z) := -r^-P n {z) exp(in^). 

^ irt'Tr I 


ITITT _ 

n£Z 

With this definition we have 

(132) t/[? 7 ](a;, 0 ) = t\(x)e 2 and J7”[ 77 ](a:, 1) = 0 for all a: £ (0,L), 

(133) &L7[i)](a;,0) = 0, 

(134) di vU[fj] = 0 on S2i. 

Consequenlty thanks to (11331) . we deduce classically that Uh [ 77 ] defined by (I129D belongs to H 2 (Jlh). 
Moreover, for all (j, l) £ N 2 s.t. max(j, l) < k + 1, we have, because of the x-orthogonality of the 

. , TTX . 

basis x i->- exp(in—), that 

Lj 

||^T;L 2 (( 0 ,L)x( 0 , 1))|| 2 < cY, [ n n\z)\ 2 dz 


n£Z 1 


< c Ei 


l 2( j+0- 3 |jj n i 2 


\Q (l \z)\ 2 dz 


V.nEZ 
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Consequently, when 77 £ Ha 1//2 (0 ,L), we get that U[f]\ £ H*{p,\) and that the following estimate 
is satisfied 

(135) < ||Q;U fc + 1 (0,l)||||77;^ fc -^(0,L)||. 

Up to the change of variable which depends only on A and is thus bounded by a constant which 
depends only on R 0l we obtain that U\ also satisfies continuity estimates similar to (1126I) - (U27I) - 
(112811 . This ends the proof. □ 

Remark 6. An alternative construction for Uh reads: 

= *///■) 

where Bh is defined by With this construction, we may exploit the better continuity estimates 

for 77 1 —y Uh[rj\ (gain of 1/2 derivative, see (|135l ) ). Yet, this construction requires more regularity 
of h than the construction given in Propostion 1 1 61 

B.2. Regularity identities. 

Lemma 17. Let h £ ff 2 (0,T;L 2 (0,L)nI 2 (0,T,]l 4 (0,I)), such that h~ x £ L°°((0,T) x (0,L)). 
For any triplet ( w,q,b ) satisfying the regularity assumptions of Definition^ namely 

(136) b £ H\ 0, T; L 2 ( 0, L)) n L 2 { 0, T; ff # 4 (0, L)), 


(137) w£Hl(Q T ), X7 2 w£L 2 /Q t ), q £ L 2 (Q T ), X7q£L 2 /Q T ), 

and such that w(x,h(x,t),t) = dtb{x,t) £ L 2 0 ( 0 ,L) and divw = 0 , the following identities are 
satisfied: 


t r-L 


(138) / divojro, q) ■ (d t w + A • Vw — w ■ VA) = — 

JQt 


0 Jo 


4>{w, q, h)d tt b 

+ /i / |Vw| 2 (f) - /i / |Vtc| 2 (0) + [ D(w ) : ([VA] T Vu; + VA[Vw;] T - D(w ■ VA)) , 

JT(f) Jt( 0) JQ t v ' 


and 


n L ft fh ft fL 

d xxxx b d tt b - (3 / / d xx bd tt b~ 7 / / d xxt bd tt b 
Jo Jo Jo Jo 


-\d tx b(x,t)\ 2 dx - pd t b(x,t)d xx b(x,t)dx - ad tx b(x,t)d xxx b(x,t)dx 


1 0 


P f f \d tx b\ 2 - a f [ L \d txx b\ 
Jo Jo Jo Jo 


-j\d tx b(x, 0 )| 2 dx + Pd t b(x, 0 )d xx b(x, 0 )dx - ad tx b(x, 0)d xxx b(x, 0 )dz. 


Proof. The second identity, involving b only, is straightforward since applying classical interpo¬ 
lation arguments (see in particular [JO] Theorem 3.1, p. 19]) the regularity (11361) of b implies 
that 

(140) b £ C°( 0, T; R 3 (0, L)) n H 1 (0, T; H 2 { 0, L)) n C 1 (0, T; H}( 0, L)) 

We now explain how one can obtain the first identity. Assuming that (w, b) are regular enough 
we have 


divcr(uj, q ) • (dtw + A • Vw — w • VA) = 


// 

JO Jy 


0 Jy—h(x,s ) 


cr(w, q)n-(dtw + A-X7w — w-X7A)dldt — 


0 J F{s) 


cr{w,q) : S7{dtw + A- Via — w- VA). 
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But by the assumption satisfied by (w, b) on the boundary y = h(x, t) and by construction of A 

on y = h(x , t). 


w-VA = W2d2A = 0 , 

d t w + A • Vw — w ■ VA = d u b, 

Combining this with the definition Q of <fi, we have 

r t r rt 

ct(w, q)n ■ (dtw + A • Vw — w ■ X7A)dldt = — 


0 J y=h(x,s ) 


</>(w,q,h)d tt b • 


o Jo 


On the other hand, we have by construction since divw = div A = 0 

div ^ d t w + A • Vw — w ■ vAj = 0 . 

Consequently, recalling that T(t) moves along the characteristics of A, we obtain the following 
equality 


'Qt 


cr(w, q) : V(d t w + A • Vw — w ■ VA) = 2 fi D(w) : D(d t w + A • Vw — w ■ VA)) 


= /i [ \Dw\ 2 ~ /x [ \Dw\ 2 (0) + 2fi [ D{w) : ( [VA] T Vw + VA[Vw] T — D{w ■ VA)') . 

Jr(t) Jt{ o) JQ t v 1 

As already noted, thanks to the transverse motion on the fluid domain, to the fact that the trace 
of w on the moving boundary is colinear to e 2 and to the divergence free property of w, we have 



Consequently, we deduce that (11381) holds true for any regular enough couple (w,b). Even 
though all the terms in the final identity are well-defined for vector-fields w satisfying (11371) (we 
already noticed in Remark [T| that with the regularity (11371) we have that t i-a |Vw| 2 (-,t) is 
continuous) we need here to get more precise on the meaning of boundary terms in integration 
by parts. Note in particular that, with (11371) . we only have dtw € L 2 (Q t ) so that its trace on 
y = h(x,t) is not well-defined (though it is divergence-free). 


There are several ways to overcome this difficulty. For instance, one may note that the compu¬ 
tations above are completely rigorous if we assume further that V<9 t w £ L 2 (Q t ). We may then end 
the proof by a density argument. Indeed, reproducing computations in the proof of Proposition [2] 
we have that the regularity (1136D - (I140L also satisfied by h, ensures that the transformation an-Mi 
defined by: 

»< W ) = (w) e a. 

realizes an homeomorphism between the set of w with regularity (11371) and the set: 

U:=H]{Sl r x (0, T)) n T 2 ((0, T); H 2 ^)). 

Furthermore this mapping exchanges the subset of w and w whose (space)-gradients are square 
integrable on the time/space domain. As U C C([0, T]-H^(Cl i)), we note also that if w n converges 
to w in U then t >->• Jjruj |Vw n (-,t)| 2 converges to t >-> fjrm |Vw(-,t)| 2 in C([0,T]). Consequently, 
we can pass to the limit in the identity (U38|l for a sequence of w n such that the associated w n 
converges in U to w. 


We can then decompose any w £ IA into 

w = wo + U[b\ ( with U defined in the proof of Proposition [TB] ) 

and approximate b by considering its (L 2 -)orthogonal projection on the first eigenmodes of the 
beam operator ( d xxxx ) and wo by considering its (L 2 -)orthogonal projection on the first eigenmodes 
of the Stokes operator with homogeneous Dirichlet boundary conditions on z = 0 and z = 1. 


□ 
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